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ABSTRACT 



When investigating the detection performance of a passive 
homing torpedo used against shallow draft surface ships, 
certain environmental factors such as the rough sea surface 
and the bubble-dominated inhomogeneous layer near the sea 
surface have to be considered. This thesis attempts to gain 
some insight into the behavior of a homing torpedo system 
during its critical attack phase, as well as getting some 
indications of the relative importance of the scattering 
mechanisms and the induced tactical limitations. An ideal- 
ized propagation model was used as reference of comparison. 
For a given sea state and target speed the results stress 
the importance of low operating frequency as well as a high 
maximum turn rate. They also point to the importance of 
having a search depth below the bubble-dominated subsurface 
layer, and a variable speed capability during the torpedo's 
attack phase. 



4 



TABLE OF CONTENTS 



I . INTRODUCTION 12 

II. SCENARIO AND TORPEDO RUN GEOMETRY 15 

III. OCEANOGRAPHIC BACKGROUND MATERIAL 17 

IV. THE PASSIVE SONAR EQUATION 21 

V. REFERENCE MODEL 22 

A. INTRODUCTION 22 

B. IDEALIZED TRANSMISSION LOSS MODEL 22 

C. NOISE SOURCE MODEL 23 

1. General Characteristics of Noise Sources — 23 

2. The Noise Source Model 26 

D. PASSIVE MODE RECEIVER CHARACTERISTICS 31 

1. Assumptions 31 

2. Derivations 32 

E. REFERENCE DETECTION PERFORMANCE 38 

1. Introduction 38 

2. 60 kHz Case 38 

3. 30 kHz Case 40 

VI. THE EFFECT OF SURFACE SCATTERING 42 

A. OCEANOGRAPHIC DESCRIPTION OF THE SEA SURFACE — 42 

B. SCATTERING THEORY 45 

C. GEOMETRICAL SHADOWING 57 

D. ESTIMATING THE SURFACE SCATTERING EFFECT 61 

E. CONCLUSION AND DISCUSSION 67 



5 



VII. THE EFFECT OF SCATTERING AND ABSORPTION FROM THE 

SUBSURFACE OCEAN LAYER 69 

A. GENERAL SCATTERING THEORY 69 

B. ABSORPTION MODEL 70 

C. THE COHERENT INTENSITY CASE 75 

D. THE INCOHERENT INTENSITY CASE 82 

E. SUMMARY AND DISCUSSION OF THE BUBBLE 

ATTENUATION 83 

F. THE REFRACTION BY BUBBLES 85 

VIII. THE TURN RATE LIMITATION 89 

IX. CONCLUSIONS AND RECOMMENDATIONS 99 

APPENDIX A: Detailed Oceanographic Background Material-102 

APPENDIX B: Surface Scattering TI 59 Program 110 

APPENDIX C: Bubble Dynamics TI 59 Program 116 

APPENDIX D: Numerical Integration TI 59 Program 121 

APPENDIX E: Turn Rate Limitation TI 59 Programs 124 

APPENDIX F: Figures 132 

LIST OF REFERENCES 174 

INITIAL DISTRIBUTION LIST 176 



6 



LIST OF TABLES 



I. Relationship between Sea State (SS) , Wind Speed, and 

Expected Significant Wave Heights (H^) 18 

II. Source Level Data and Calculation for a Blade Surface 

Cavitating Propeller 29 

III. Detection Parameters and Ranges for f=60 kHz 86 

IV. Detection Parameters and Ranges for f=30 kHz 87 

V. Frequency Distribution of VJind Speed in Percent per 

Year at Weather Stations Along the Norwegian Coast 104 

VI. Frequency Distribution of Significant Wave Height in 

Percent at the Weather Station "Andenes" 106 



) 

7 



LIST OF FIGURES 



1. Torpedo Search and Attack Geometry 132 

2. Location of Weather Stations along the Norwegian 

Coast 133 

3. Average Monthly Wind Speed in Beaufort from the 

Weather Station "Andenes" 134 

4. Histogram of Significant Wave in Percent per Year 

from the Weather Station "Andenes" -135 

5. Resonant Bubble Density at 12 kHz as a Function of 

Depth 136 

6. Resonant Bubble Density at 38 kHz as a Function of 

Depth 137 

7. Resonant Bubble Density at 120 kHz as a Function of 

Depth 138 

8. Ambient Noise Level Curves 139 

9. Typical Sound Speed Profiles in Norwegian Coastal 

Waters 140 

10. Worst Case Ray Path During Winter 141 

11. Worst Case Ray Path During Summer 142 

12 . Frequencies of Wind Directions in Percent for 

Stations along the Coast from "Hillesoy" to "Ona" 143 

13. Frequencies of Wind Directions in Percent for 

Stations along the Coast from "Myken" to "Furuholmen"-144 

14. Absorption Coefficient in dB/m as a Function of 

Frequency 1 4 5 

15. Average Radiated Spectrum Level for Surface Ship as 

a Function of Speed in kts 146 

16. General Noise Spectrum for a Cavitating Propeller 147 

17. Spectrxjm Level as a Function of Total Gas Content 148 

18. Measured Model and Full Scale Noise Spectra 149 



8 



m 




19. 



20 . 



21 . 



22 . 

23. 

24a. 

24b. 

25. 

26. 
27. 



28. 



29a. 

29b. 

30. 

31. 

32. 

33. 



34. 



35. 



36. 

37. 



38. 



39. 



40 . 



Square Law Detector Scheme 150 

Required Input S/N Ratio and BT Product for Various 
Operating Probabilities 151 

Probability of Detection Versus Input S/N Ratio for 
Various Thresholds 152 

-201ogR-aR Versus R for a Frequency of 60 kHz 153 

-201ogR-aR Versus R for a Frequency of 30 kHz 154 

Geometry of the Sea Surface Scattering 155 

Specular Scattering Geometry 156 

Rough Surface Shadowing Geometry 157 

Calculation Scheme of the Surface Scattering Effect 158 

AIL Versus R 159 

Geometry of the Subsurface Propagation Model 160 

Chains of Successive Scattering 161 

Scattering Path Going Through the Same Scatterer More 
Than Once 151 

Scattering Pattern 162 

a and a for 60 kHz 163 



e a 

a and a for 30 kHz 

e a 

Interpolated Bubble Data 

a^(a)n(a)da for 60 kHz 

a^(a)n(a)da for 30 kHz 

a^(a)n(a)da for 60 kHz 

a (a)n(a)da for 30 kHz 

Pursuit Homing Geometry 

I $K/V^ 1 as a Function of (j> with Parameter p 
Hit Criterion Geometry 



164 

165 

166 

167 

168 

169 

170 

171 

172 



9 



173 



41. 



AOB-Limitation Versus Ship Speed with Parameter Maximum 
Turn Rate 



10 



AC KNOWLEDGMENT 



The author wishes to gratefully express his apprecia- 
tion to his advisor. Professor Karlheinz E. Woehler, for 
his advice, encouragement and guidance in the preparation 
of this thesis. 



11 



I. INTRODUCTION 



The following analysis examines several factors that 
limit the detection performance of a passive homing torpedo 
with the mission objective of countering shallow-draft 
targets in Norwegian coastal waters. Generally, these 
factors can be divided into three main groups: 

-Environmental factors in the ocean 

-Electrical mechanical and hydrodynamical factors in 
the torpedo system 
-The users tactical situation 

The factors that are generated in the ocean itself are the 
principal subject of this analysis. As our interest is 
confined to the layer immediately below the surface of the 
ocean, later called the subsurface ocean layer, the main 
factors affecting the sound propagation are: 

-Scattering and absorption due to the bubble-dominated 
inhomogeneous subsurface ocean layer. 

-Scattering from the rough sea surface 
The concentration of air bubbles and the roughness of the 
sea surface are determined by the windspeed. The effective- 
ness of both these scattering mechanisms depends on the 
frequency of the incident wave and the geometry of the 
source and receiver. The following analysis is limited to 
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high frequencies in the region of 30-60 kHz which are 
characteristic of existing torpedo systems. At high fre- 
quency and low grazing angles for the incident and 
received signals, the phenomena of "shadowing" of the 
surface by other parts of the boundary occurs. Under 
these conditions, the effect from the inhomogeneous sub- 
surface ocean layer becomes increasingly important. Ob- 
viously both the above mentioned scattering mechanisms are 
present simultaneously. Often these two effects cannot be 
resolved either theoretically or experimentally, as any 
signal with a finite duration will be scattered from the 
space near the sixrface simultaneously with that from the 
sea surface itself. 

In order to adequately describe the scattering mech- 
anisms, this analysis starts with a presentation of the 
oceanographic background material for : 

-Typical windspeed and wave height 

-Typical sound speed profiles 

-Density and distribution of air bubbles in the sub- 
surface layer 

-Statistical description of the sea surface 
The analysis proceeds by separately estimating the effect 
of : 

-Scattering from a randomly rough surface 

-Scattering and absorption caused by an inhomogeneous 
sxibsurface layer. 
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and comparing their relative importance. The method 
employed for these estimations is an approximation that 
is a combination of both ray and wave theories. Ray 
methods are used to follow the acoustic signal from its 
source to the vicinity of the scatterer. Wave theory is 
used to calculate the actual scattering process. Finally, 
ray theory is again used to follow the scattered signal 
to the receiver. 

An idealized propagation model consisting of an iso- 
tropic stratified medium will be used as reference of 
comparison. This model is founded on: 

-A noise source from a cavitating propeller blade. 

-The operational characteristics for a square law 
detector (ROC-curves) . 

-A transmission loss model based on geometrical 
spreading and absorption losses in homogeneous sea 
water. 

In this analysis, the passive sonar equation is used to 
predict the performance of the homing system. The detec- 
tion range encountering the two scattering effects will 
be obtained from the sonar equation and compared to the 
detection range based on the reference model. Thus, the 
difference in ranges at which the homing device just 
acquires the target with and without scattering is a 
measure of efficiency. 
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II. SCENARIO AND TORPEDO RUN GEOMETRY 



In the Noirwegian coastal waters, the primary mission of 
a torpedo system is to counter an amphibious force consist- 
ing of escorts, supply ships, and shallow-draft landing 
crafts. Typical characteristics of these three ship 
types are as follows: 

Supply ships; 

-Displacement 5000 tons 
-Length 100 m. 

-Draft 6 m. 

-Speed 15-20 kts. 

Escort ships: 

-Displacement 2000-3000 tons 
-Length 85 m. 

-Draft 3 m. 

-Speed 35 kts. 

Landing craft; 

-Displacement 1000 tons 
-Length 80 m. 

-Draft 2 m. 

-Speed 18 kts. 

In order to simplify this analysis, moderate sea states (SS 3) 
are assumed. Since "moderate" wave heights of two meters 
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are appreciable when compared with the two-meters draft 
of the landing craft, the possibility of an acoustic torpedo 
impacting the target at a depth of two meters is very remote 
without the use of an influence exploder. 

Two relevant search and attack schemes will be con- 
sidered. These are illustrated in Fig. 1 together with the 
operation of the influence exploder. From target valida- 
tion to completion of terminal attack, the torpedo contin- 
uously tracks in the azimuth plane. In Case A, ascent is 
inhibited after enable. For Case B, ascent is inhibited 
after the torpedo reaches terminal attack depth. 

An assumed attack depth of six meters is consistent 
with the activation range of influence exploders and is 
deep enough to preclude wave or "free surface" induced 
disturbances of the torpedo. Success of the attack depends 
primarily on the availability of maintaining azimuth-plane 
steering to within a short horizontal range of the target, 
and the subsequent operation of the influence exploder. 

Case A is of particular interest- to this analysis, as 
both the search and attack-depth are within the subsurface 
layer . 
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Ill- OCEANOGRAPHIC BACKGROUND MATERIAL 



The oceanographic background for predicting typical 
and extreme conditions of 
-wind speed 
-wave heights 

-bubble densities and distributions 
-ambient noise versus self noise 
-sound speed profiles 

are outlined in detail in Appendix A. Even though most 
data have beneral validity for Norwegian coastal waters, 
the region above 68°N are of particular interest. Thus, 
a typical area combining open as well as confined waters 
can be represented by "Andfjord" at 70°N, where the oceano- 
graphic conditions can be related to the weather station 
"Andenes," see Fig. 2. 

Figures 3 and 4 [Ref. 1] show average windspeed and 
the occurrence of significant wave heights as a function 
of time of year at weather station "Andenes," respectively. 
The bulk of data is centered around a windspeed of 
Beaufort:4-5 (11-21 kts) and SS:3-4 (significant wave heights 
1-2 m) . Table I gives the relationship between SS, wind- 
speed and expected significant wave heights. 
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TABLE I 



RELATIONSHIP BETWEEN SEA STATE (SS) , WIND SPEED, 
AND EXPECTED SIGNIFICANT WAVE HEIGHTS (H ) 

S 



Class 

Number 



Significant Wave Height 
in m 



Beaufort 

Scale 



SS 



H 



s 



0 


0 


0 


1 


0 - 0.1 


1 


2 


0.1 - 0.5 


2 


3 


0.5 - 1.25 


3-4 


4 


1.25 - 2.5 


5 


5 


2.5 - 4.0 


6-7 


6 


4.0 - 6.0 


8 


7 


6.0 - 9.0 


9-10 


8 


9.0 - 14.0 


11 


9 


>14 


12 
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This, together with the low probability of having an 
amnibious operation occurring in high sea states (SS>5) 
justifies the assumption of moderate sea state with wind- 
speed in the region of 12 kts and wave heights of 2 m. 

The bubble data distribution taken from Ref. 2 was 
obtained in the area "Tromso" - "Bjorndya" during the 
period June-November 1978. These data correlate very well 
with a larger body of data obtained by H. Medwin [Ref. 3]. 

Figures 5, 6 and 7 [Ref. 2] show the density of resonant 
bubbles as a function of depth with windspeed as parameter 
for the 12,38 and 120 kHz. As seen from these figures, 
the number of resonant bubbles are an increasing function 
of frequency and windspeed, and a decreasing function of 
depth. Below a depth of approximately 15 m, the number of 
bubbles is negligible for the windspeed of interest. 

The effect of SS (windforce) on the ambient noise level 
is given in Fig. 8 [Ref. 4] . Shallow coastal Norwegian 
waters are typically 5-10 dB noisier than the corresponding 
deep water. However, great variability caused by local 
ship traffic, fishing fleet activity, marine life and 
local wind conditions makes ambient noise prediction diffi- 
cult in these areas. This means that for accurate modeling, 
ambient noise prediction have to be done at each location 
as its level is both site and time dependent. However, 

Fig. 8 shows that for frequencies higher than 50 kHz, the 
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effect of windforce on the ambient noise level decreases 
to a lower bound determined by the thermal agitation. 

Based on the above discussion and experience related to 
noise levels for torpedo systems, the self noise will be 
assumed to be dominant through this analysis. 

Figure 9, obtained from Ref. 5, shows that the sound 
speed profiles usually encountered in the area of interest 
results in extremely difficult sonar conditions. This is 
illustrated in Figs. 10 and 11, which show worst-case ray 
path derived from Fig. 9. In addition, the presence of 
bubbles in the subsurface layer causes the sound speed to 
be a function of frequency. The above two factors may 
frequently be the ones limiting the detection range of the 
torpedo. These effects can be minimized by selecting an 
appropriate search depth for a particular sound speed 
profile. In addition, for Case B the corresponding curved 
homing trajectories in the pitch plane give an error in 
apparent range to target. This effect will normally be 
taken into account by devising appropriate attack logic 
which is outside the scope of the present analysis. 
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IV . THE PASSIVE SONAR EQUATION 



A measure of efficiency for a passive homing torpedo 
is the detection range obtained by solving the passive 
sonar equation for broadband noise: 

SL+lOlogB-TL (geom) -aR+DI-NL-DT=0 (1) 

where 

SL = spectral level of the broadband noise 

radiated by the target (in dB re luPa/lHz at 
Im) . 

R = detection range (in ra) . 

a = attenuation coefficient at the center frequency 
(in dB/m) . 

DI = receiving sensitivity (directivity index) 

(in dB re luPa) . 

NL = noise level at the receiver in the bandwidth B 
(in dB re lyPa) . 

DT = detection threshold; the signal to noise ratio 
at the transducer output required for a de- 
tection probability of p^^ and associated false 
alarm probability Pp^ (in dB) . 
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V . REFERENCE MODEL 



A. INTRODUCTION 

In order to produce the reference for the analysis 
the sonar equation is solved assuming ideal free-field 
conditions, a simple noise source model, and a generalized 
square-law detector. 

B. IDEALIZER TRANSMISSION LOSS MODEL 

Because the presence of refraction, scattering, and of 
ocean boundaries, free-field conditions associated with 
homogeneous (isovelocity) and unbounded medium seldom 
exist in the sea. However, as a basis for comparison, 
the ubigutuous spherical spreading law plus an added loss 
term due to "normal absorption" can be used as a reference 
model for measuring the effects of the previous mentioned 
scattering and absorption mechanisms. Thus, the reference 
transmission loss model can be expressed as; 

TL=201ogR + aR (2) 

where the absorption coefficient, expressed in dB/m, can 
be obtained from Fig. 14 taken from Ref. 8. 
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C. NOISE SOURCE MODEL 



1 . General Characteristics of Noise Sources 

Soxand is generated in a fluid medium by any process 
that causes an unsteady pressure field. Physically processes 
that can cause an unsteady pressure field include: 

-Pulsation of a boundary surface of the medium 
-The action of a nonsteady source on the fluid 
-Turbulent motion in the fluid 
-Oscillatory temperatures 

It can be shown, e.g.. Ref. 10, that each source mechanism 
mathematically corresponds to a dominant order of multipole. 
If all sources are of such a nature that their time variation 
can be described by a Fourier Integral, it can be shown 
[Ref. 10] that the Helmholtz Equation is 

^ -4-irf^(x) (3) 



9Q,,(x) 

OJ 

3t 



+ 



V*F (x) - 
0) 



Term 1 Term 2 

where the right hand side describes distributed 
The terms on the right hand side of Eq. 
following interpretations: 

Term 1: mass injection 

Term 2 : external force 

Term 3: turbulent shear stress 



3x^3Xj 
Term 3 

source terms. 
(3) have the 



) 
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In the long distance and long wavelength approximation, 
it can be shown that the mass injection term gives rise to 
a simple source; a zero order pole called a monopole. The 
monopole radiates omnidirectional with no angular dependence. 
At large distances the pressure field from the monopole 
radition is that of a point source. Examples of this are: 
-Pulsating bubbles 
-Cavitation 

The external force, in the long distance and long wave 
length approximation, is associated with a dominant dipole 
which has a cosine directional pattern. Examples of this 
type of radiation is that caused by the vibratory motion 
of an unbaffled rigid body. 

Radiation from turbulent shear stresses is charac- 
terized by a lowest order term of quadrupole nature. 

The efficiency of the source terms decreases with 
increasing dependence on the spatial derivatives. This 
can be understood when recognizing that wave functions of 
the general form f(x-ct) have a time derivative 

|-^f(x-ct)| = cf^(x-ct) (4) 

which is magnitude c (sound speed) greater than the spatial 
derivative 

l^f(x-ct)! = f^(x-ct) (5) 
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other factors being equal, the radiation from an external 
force is small compared to that from mass injection, and 
that from turbulent shear stress is the smallest; therefore, 
moncpole radiation is the dominant term. 

Propeller cavitation, when it occurs, is usually 
the dominant noise source for any marine vessel. Submarine 
and torpedoes often operate at a depth great enough to avoid 
cavitation. Surface ships, on the other hand, generally 
have well developed propeller cavitation with the result 
that their radiated spectrum from 5 Hz to 100 kHz is 

controlled by this source. 

The basic phenomena of cavitation combined with 

propeller hydrodynamics give the fundamental characteristics 
of propeller cavitation noise. An excellent qualitatively 
discussion of this can be found in Ref. 9:Chs. 7 and 8, 
from which the following is extracted: 

Propeller blades are rotating twisted wings 
that produce hydrodynamic forces. Depending on operating 
conditions, they experience cavitation on a number of 
different places. Of these there are three prominent types: 
-Tip vortex cavitation 
-Hub vortex cavitation 
-Surface blade cavitation 
In addition to the two types of vortex cavitation, 
there normally are two types of blade surface cavitation: 
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-Back: driving face 

-Front: suction surface 



Of all kinds of propeller cavitation, surface blade cavi- 
tation on the suction surface is normally the most noisy, 
while hub vortex cavitation is the least noisy. 

2 . The Noise Source Model 

Due to lack of recorded and available noise data 
from the target in question, the noise source has been 
generalized on the basis of the following discussion and 
assumptions : 



blade cavitating propeller operating in a good to poor 
wake; surface cavitation will be assumed to be dominant. 



theory for cavitation noise, where dimensional analysis 
is combined with the basic results of cavitation theory 
that the acoustic pressure is proportional to the product 
of the collapse pressure of the cavities and the volume 
of cavitation produced per unit time. From this synthesis 
it is found that the total acoustic intensity varies as: 



The noise source will be build up around a surface 



D. Ross [Ref. 9] has developed an approximate 



I = K 



ti 




( 6 ) 



where 



r = distance of the hydrophone from the source 



K. . = the cavitation inception parameter 
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= blade tip speed 

~ t)lade tip speed for inception of cavitation 

This expression shows that propeller cavitation noise power 
is proportional to the total number of blades, b, the blade 
chord, s, and to the propeller diameter, D, and is a function 
of the tip speed with the dependence on the tip speed being 
the strongest. The different blade surface sections where 
cavitation exist are uncorrelated and the radiated noise 
is treated as a single monopole radiation so that at a 
distance r>>a (where a is the characteristic dimension on 
the source region) the radiation is similar to that of a 
point source with no angular dependence. 

Submarines and torpedoes with centerline propellers 
have a relatively symmetric inflow condition. Surface ship 
propellers, in contrast, operate under nonxiniform inflow 
conditions. Circumferential wake variation causes the 
radiated sound to be amplitude modulated at the blade rate 
frequency. Furthermore, slight physical difference 
between the blades produces modulation at the shaft rate 
frequency. These amplitude variations gives a very 
distinct characteristic to the radiated noise that can be 
used for classification purposes to reduce the probability 
of false alarm. 

The most complete source of data on surface ship 
radiated noise are measurements made during WWII, reported 
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in a compendium issued by the U.S. Office of Scientific 
Research and Development (OSRD) in 1945 and declassified 
in 1960. When these data are examined the radiated noise 
is found to depend on tip speed and the number of propeller 
blades with little dependence on the other variables. For 
surface ships near cruise speed, the source level for 
frequency over 100 Hz can be written as: 

SL=SL'+20-201ogf ; f>100 Hz (7) 

where 

f = frequency in Hz 

SL' = overall source level in dB re lyPa. 

The overall level can be expressed as: 

SL'=175+601ogU /25 +101ogb/4 (8) 

where 

U = TTnD (9) 

n = rotational speed (rpm) 

D = diameter of the propeller (m) . 

The above expressions are used as the basis for the noise 
model with the following input data: 

n = 300 rpm for maximum cruise speed of 15 kts. 
n = 180 rpm for a cruise speed of 10 kts. 

D = 2 m. 
b = 5. 

The resulting noise spectrum, in dB re lyPa at 1 m, 
as a function of speed in kts. are tabulated in Table I*- 
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SOURCE LEVEL DATA AND CALCULATION 
A BLADE SURFACE CAVITATING PROPELLER 
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and plotted in Fig. 15. A one sigma region (5 dB of un- 
certainty) is incorporated in Fig. 15. 

The above noise spectrum estimation agrees very 
well in the high frequency limit, with more recent studies 
by A. Ldvik [Refs. 11 and 12] . Here it is found that the 
cavitation spectra, both theoretically and experimentally, 
can be divided into four frequency regions, as illustrated 
in Fig. 16. 

Region I is dominated by noise at the blade fre- 
quency and its harmonics. The emitted sound is caused by 
the volume variation of the main cavity. 

Region II starts at the bubble frequency, which is 
the reciprocal of the lifetime of the main cavity. The 

mean power level is found to decrease with increasing fre- 
.-2.5 

quency as f 

Region III is an intermediate region. 

Region IV associated with the shock waves starts 
at the mean collapse frequency f^, given by the mean 

collapse time. The power level is found to decrease as 

-2 

f , as in Eq. 5. 

The number of gas bubbles in the water have a 
pronounced effect on the high frequency cavitation noise 
from the propeller. This is illustrated in Fig. 17, 
obtained from Ref. 12, where the power is found to decrease 
as much as 40 dB with increasing gas content. 
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Scaling laws are developed [Ref. 12] for each region 
based on a series of models and full scale measurements. 

These laws depend on the dynamic pressure induced by the 
propeller, the model ratio, and the gas content of the water. 
The full scale measurements were performed in cooperation 
with the Royal Norwegian Navy and the Marine Institute of 
Norway. The model experiments were performed in the largest 
cavitation tunnel at the Ship Research Institute of Norway. 

In summary, the scaling of cavitation noise was 
demonstrated to be a useful tool in predicting a full scale 
cavitation noise as shown in Fig. 18 [Ref. 12] which compares 
measured noise spectra for the model and full scale 
measurements . 

For the high frequency region, the source levels 
are of the same magnitude as predicted by the WWII empirical 
formula. 

D. PASSIVE MODE RECEIVER CHARACTERISTICS 
1. Assumptions 

For receiver characteristics assume a square law 
detector with a center frequency f-60 or 30 kHz and a band- 
width B. The detection scheme is shown in Fig. 16. The 
principal assumptions employed in the derivation are as 
follows : 

-Gaussian signals in Gaussian noise 

-Frequency independent signal and noise spectra 



31 



-Integration time T is sufficiently long to 
permit application of the central limit theorem. 



2 . Derivations 

The detector input r(t) is assumed to be a zero- 
mean Gaussian process composed of noise alone or signal 
plus noise expressed by the two well known hypotheses 

H^r(t) = n(t) (10) 

Hi r (t) = s (t) +n (t) 

where : 

n(t): noise signal 

s (t) : signal . 

The two signals s (t) and n(t) are assumed to be independent. 

Assume that the spectral shape of s(t) and n(t) 
are the same, such that and only differ in the total 




Schematic of detector-smoother 



and 

X (t) = r^ (t) (11) 
Furthermore, let the noise variance be normalized to unity 
(for convenience) and the signal variance by denoted by 



Var [n (t) ]=1 
Var [s (t) ]=a^ 
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( 12 ) 



Because of the assumed similarity in the spectral shapes, 
the autocorrelation functions are 

J^n ] “P N(f) is the noise power (13) 

spectral density. 

Rg ( T ) = p ( X ) a ^ 

Rn(T)= p(x) 

(l+a^)p(T) 

Furthermore, assume that the integration time T is long 

enough so the central limit theorem holds, implying that 

K also is a Gaussian random variable. 

This yields that the probability density function 

of the output variable and hence the detection and false 

alarm probabilities are completely determined once the 

mean and the variance of K are derived. 

If a process V(t) is wide-sense stationary, then 
T T 

E[V{t)]=E^ yv(t)dt yE[V(t)]=v (constant) (14) 
0 0 



Thus, assiaming that r(t) is a wide-sense stationary process. 

E[K]=E[x(t) ]=E[r (t) ]=l+o^ (15) 

and similarly 



TT 

Var[V]=E[V^]-{E[V] }^= A- Jj' E[V(t) V(s) ]dtds-v^ (16) 

00 



TT TT 

Var[V]=-i- yy^ [R^(t-s) -v^]dtds=~ // Cov^(t-s)dtds 
00 00 
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where Cov^(t-s) is a covariance function. 
Then letting 



X = t-s 

^ = t+s, and substituting into Eq. (7) yields 

Var[V]=i^ov^(T)4^ 

o 

T 

Var[V}=^ J [l~-^-^]Cov^ (t) dx 
-T 




Consequently 



Var[K]=^ 



-Ml 



[l-^^]Cov^(x)dx 



-T 



Then, evaluating the covariance function from the auto 
correlation function 

(x)=E[x(t) x(t-x) ] = E [r^ (t) r^ (t-x) ] 

Since r(t) is Gaussian, the above fourth moment can be 
expressed as product and sums of second moments; 



R (x)=r^(t) • r^ (t+x) + 2r(t)r(t-x) • r(t)r(t-x) 
=rJ(0) + 2r2(t) 

R^(x) = (1+a^) ^+2 (1+a^)^ p ^ (x) 

Thus, the covariance function is 
Cov^(x)=R^(x)-{E[x(t ) ] V 

= (l+a^) 2+2(l+aM pMx)-(l+a2) 2 



(17) 



(18) 



(19) 



( 20 ) 
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Cov^(t)= 2 (1+a^) ( t) 

Inserting Eq. (12) into Eq. (9) yields 

T 

Var[K]=-i- J {[1-1^1 2(1+o^)^o2(t)) dt 

-T 

Var[K]=iii^m j [1-M] 

-T 

If T is large compared to the correlation time 
TB>>1 

then we can substitute the limit for Eq. (22) by 



Var [K]= 



2 ( 1 + 02)2 
T 



00 



— 00 



(t) dt 



(23) 



2 (l+a2) 2 
T 



00 



— 00 



N2 (f)df 



If we further make the assumption that the signal and noise 
have ideal flat bandpass spectra: 

N(f)= f 1/2B, f-B/2<f<f+b/2 (24) 

[o, otherwise, 

Inserting Eq. (24) in Eq. (23) yields 



Var [K]= 



(1+g^) " 
BT 



(25) 



The probability density function for the output of the 
detector have the following form 
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) 




Here w is the outcome of all possible signals. The false 
alarm probability is obtained by integrating the conditional 
probability p(w|H^) over the outcome space for which to 
choose 



00 




Similarity, the miss probability=l-detection probability 
is obtained by integrating the conditional probability 
p(w|Hj^) over the outcome space for which to choose 

P„ = l-Pj, = / p{w|H^)dw = Q(-^) (27) 

.00 

Further defining the input and output signal-to-noise 
ratios as 

2 

S/N (input) = ^ (28) 



(Pi-Mo) ^ 

S/N (output) = (29) 

0i 
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For the square law detector where 

= E[n^(t)] =1, yj = E[r^(t)] =E[K] =l+a^ (30) 

Og" = Var[n^(t)] = = Var[r2(t)] 

= Var[K] = ' 

L J g,p 

The output signal-to-noise ratio: 

S/N (output) = — ( 31 ) 

r (l+oM (l+a^)2 

L Qrj, J 

The input signal-to-noise ratio: 

S/N(input) = (32) 

The probability of false alarm: 

Pp^ = Q[/BT(j-l) ] (33) 

The probability of detection: 

p = 1-Q[/BT (- — a^ ) ] (34) 

1+a 

The Equations (32) , (33) , and (34) are plotted for a 

variety of realistically encountered conditions. Figure 17 
gives BT versus S/N (in) (Identical to the detection 
threshold DT) for various combinations of Pj^ and 

Figure 18 gives p^^ versus S/N (in) for various 
combinations of the threshold j and Pp^* This constitutes 
the ROC-curves for the generalized square law detector. 
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E. REFERENCE DETECTION PERFORMANCE 

1 . Introduction 

As both the scattering and absorption are frequency 
dependent, it is necessary to have reference models for both 
60 and 30 kHz. 

2 . 60 kHz Case 

The range dependent portion of the sonar Equation (1) 

-201ogr-aR 

is plotted in Fig. 19 for a=0.021 dB/m taken from Fig. 14. 

Figures 15, 17 and 19 are then used to estimate the 
detection range: 

(a) Select p^, Pp^^/ B, and the integration time T. 
The detection threshold DT=S/N (input) is then 
found from Fig. 17. 

(b) Select the speed of the target and find the SL 
from the noise source model (Fig. 15) . Then, 
reasonable values for the receiver sensitivity 
DI and the self noise level NL yields the left 
hand side of the sonar equation (1) except for 
range dependent term. 

(c) Use Fig. 19 to solve the passive equation 
for R. 

A realistic example may illustrate the above procedure. 
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(a) Entering Fig. 17 with 
Pd = 0.5 

PPA = 

B = 4500 Hz 
T = 100 msec 

yields 

DT = -6.5 dB. 

(b) Selecting a target speed of 12 kts. gives 
(from Fig. 15) SL=100 dB. Selecting a 
typical transducer sensitivity DI = -180 dB. 
Assxaming the NL to be dominated by self noise 
of typical value NL=-124 dB. This yields 

SL+101ogB+DI-NL-DT=87 dB. 

(c) Figure 19 with 

-201ogR-0. 02R=-87 
yields 

R=1200 m. for a=0.02 dB/m and f-50 kHz. 

The influence of different design parameters like self 
noise and detection threshold on the passive detection 
performance is now easily investigated by the above 
procedure. 

Although outside the main scope of this analysis, 
the above statement can be confirmed with an example. As 
seen from Fig. 18, a probability of false alarm Pp^ = 10 ^ 
implies a threshold setting = 0.9 dB above the noise level. 
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Due to the variability of the noise level together with 
the practical difficulty in accurately setting the threshold, 
a more realistic goal for the threshold would typically be 
3 dB. Going into Fig. 17 shows that the corresponding value 
for DT for p^^ = 0.5 is DT = 0 dB, with a corresponding low 

value for the p„. . 

*^FA 

Letting DT=0 dB and keeping the previous assumed 
values of SL, DI, and NL yields; 

SL+101ogB+DI-NL-DT=80. 5 dB. 

The corresponding detection range is; 

R = 1000 m, for a = 0.02 dB/m and f = 60 kHz. 

Thus, this change in threshold setting caused a decrease 
in detection range from 1200 m to 1000 m in return of a 
significant decrease in the false alarm probability. 

3 . 30 kHz Case 

In order to estimate the reference detection range 
for an operating frequency of 30 kHz, we utilize the sonar 
equation (1) . 

Assume that the receiver has the same generalized 
passive detector characteristics as in the 60 kHz case; 

DI = -180 dB 

NL = -124 dB. 

_ g 

DT = -6.5 dB, based on Pj^ = 0.5 and p^^ = 10 
However, the empirical equation (7) for the source level; 

SL=SL'+20-201ogf 
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I ^ 




shows that the source level falls off as f“^. If the 
dynamical and dimensional parameters of the propeller are 
the same, SL will increase by +6 dB when the frequency is 
reduced from 60 to 30 kHz. 

A source level of 
SL=100+5=106 dB 

gives a range dependent solution of the sonar equation 
-201ogR -aR=-87-6 = -93 dB. 

A plot of 
-201ogR-aR 

is given in Fig. 20 for an absorption coefficient a=0.01 dB/m 
taken from Fig. 14. 

Figure 20 then gives a detection range of 
R = 2400 m. 

Thus, as seen from these ideal reference calculations, 
halving the frequency gives a higher source level and a 
lower absorption loss, resulting in a doubling of the 
detection range. 
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VI . THE EFFECT OF SURFACE SCATTERING 



A. OCEANOGRAPHIC DESCRIPTION OF THE SEA SURFACE 

The roughness of the sea surface is the essence of 
the scattering mechanism. Thus, to adequately describe 
the scattering of sound from a randomly rough sea surface, 
it is necessary to formulate a suitable description of the 
sea surface from an acoustical propagation point of view. 

Generally the shape of the rough sea surface is most 
appropriately described in terms of time and spatial 
dependent random variables. However, observation of the 
ocean under the same environmental (meteorological) condi- 
tions indicates that the roughness is the same over large 
areas and for periods of at least several hours. The 
random processes responsible for the structure of the sea 
surface, therefore, can be considered stationary at least 
over periods of hours. With this assumption, the sea 
surface can be described in terms of the statistical 
description of the surface displacement function, and the 
distribution of signals reflected from the sea surface 
can then be related to this probability distribution 
function . 

The most significant statistical parameters describ- 
ing the scattering mechanism from the randomly rough 
surface are: 
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‘i| 






-the mean square slope 

-the mean square surface height 

-the correlation length. 

Optical measurements made at sea by C. Cox and W. Munk 
[Ref. 13] showed that the sea surface with an arbitrary 
wide continuous spectrum of waves is characterized by a 
Gaussian distributed surface slope. The mean square 
slope, determined from these optical measurements is 

<i;i 2 > = s'" = (3+5.12w)xl0“^ (35) 

where 

W=wind speed in m/s measured 41 ft {12.5 m) above 
the sea surface. 

The Gaussian distribution of surface slopes implies that 
the surface displacement function can be described by a 
Gaussian probability density function with zero mean 
<^> = 0 
and variance 

<;> = 0 ^ 

and Gaussian correlation function 
(t) ; (t+t) > 

The mean square height o^, is obtained by integrating the 
frequency spectrum of the fully developed sea. The fre- 
quency spectrum G is given by the commonly accepted semi- 
empirical expression of Piersom— Moskowitz [Ref. 14] as 



) 
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G{Q) 



(36) 



_ ag^ p ^Oi 4 

= exp [-6 

where 

= frequency (in s~^) 
a = 8.1 
6 = .74 

^Q= g/W (in s~^) . 

W = wind speed in m/s at 19.5m above the sea surface 
g = gravitational acceleration (in m/s^). 

This gives 



= 



G(n)dfi = 



aW 



4Bg‘ 



(37) 



For a Gaussian autocorrelation function expressed as 

’Mt) = — <; (t) ; (t+T) >= ^ ^ (38) 

where 

T = correlation length. 

The following relationship for the mean square height holds 
for sea of small roughness 



•2 _ 2o^ 



= 



, see later Eq. (60) 



(39) 



and the correlation length is thus: 
T = /2 ^ 



(40) 
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B . SCATTERING THEORY 



All real boundaries are rough for radiation with short 
enough wavelength, and the apparent roughness depends on 
the "viewing" conditions. The wave reflected by a plane 
surface has the same properties as the incident wave since 
the radiation is scattered coherently and there is a definite 
relation between the incident and scattered waves. 

A randomly rough surface, however, such as the wind 
generated ocean surface, scatters radiation in all directions, 
i.e., an illuminated area is visible from any direction. 

Heuristically there are two distinct approaches to this 
phenomena . 

1. If the boundary is rough most of the radiation is 
scattered and there is little transmission in the 
specular direction. Thus, the attenuation caused 
by the irregularities can be included in the trans- 
mission equation. 

2. If the surface is truly smooth, it can be assigned 
that the effect of the boundary is to supplement the 
original pressure field by an out-of-phase image 
contribution. For a randomly rough surface the 
reflected sound neither completely cancels the 
direct sound nor adds to give +6 dB pressure peaks 
of the interference pattern. For a rough surface, 

this supplement is a small fraction of the direct 
path. 
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The second approach will be used with the simplified 
assumption that the sea below the surface has an isotropic 
statistical description; i.e., the mean acoustic velocity 
and the mean density are assumed to be constant and have 
negligibly small mean square fluctuations. 

The estimation of the scattering is based on an approxi- 
mation method employing both ray and wave theory. Ray 
methods are used to follow the acoustic signal from the 
noise source to the vicinity of the sea surface. Then, 
wave theory is used to calculate the scattering process. 
Finally, ray theory is used to follow the scattered signal 
to the receiver. 

The geometry is given in Fig. 24 a. The origin of the 
coordinate system is at the center of the illuminated area. 
The x-y plane coincides with the mean of the rough surface 
as averaged over the illuminated area. 

The source and receiver are at distances and R2^ 
respectively, from the origin. R^ is the xz-plane and 
makes the angle 0 ^^ with the z-axis. R2 makes the angle 02 
with the z-axis and the projection of R2 on the xy-plane 
has an angle 0 ^ relative to the x-axis. 

For high frequencies R^ and R2 are much larger than 
the acoustic wavelength. Then both the incident wave 
and scattered waves can be treated as nearly plane waves. 



46 



The formulation of the scattering problem will be 
based on the Helmholtz integral which requires known 
values of the normal derivatives of the incident and 
reflected waves on the boundary. These are estimated 
by using the neuristic Kirchhoff's approximation, which 
assumes that the wave is locally reflected by a plane 
surface; i.e., an approximation restricted to a surface 
not too rough and not shadowed. 

Further, the receiver derictivity, as indicated by 
Fig. 24b, will be used to limit the surface area that is 
illuminated . 

Since this procedure is based on a detailed develop- 
ment by I. Tolstoy and C. S. Clay [Ref. 15], only the 
main points will be outlined here to bring out the assump- 
tions made and the inherent limitations of this approach. 

The development starts by considering the inhomogeneous 
wave equation of the general form: 

^ . 3^p{Xit) 

V^p(Xit) - — -4tt f(Xit) (41) 

c^ 3t^ 

where 

f(Xjt) = is a known source distribution. 

The development is based on the following initial assumptions 

-The medium is homogeneous. 

-The medium is bounded by some surface S, onto which 
an incident wave impinges. 



■) 
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-The boundary is characterized by the specific 
acoustic admittance and the shape of the boundary. 

-The incident wave is harmonic. 

The assumed harmonic source implies that f(X/t) can be 
decomposed into a Fourier integral. Furthermore, assuming 
that the solution of Eq. (41) can be decomposed in time, 
we arrive at the Helmholtz equation 



It should be noted here that L. Fortuin, in [Ref. 16] 
showed that the Helmholtz equation is not exactly correct 
for a medium with a time dependent boundary. The equation 
can, however, be used with a good approximation when the 
time derivative of the surface elevation is much smaller 
than the speed of the waves through the medixam. For under- 
water sound waves scattered by the rough sea surface, this 
means that the wind speed has to be much less than the 
sound speed; a requirement easily fulfilled for our inves- 
tigation. 

Green's method allows the solution of this linear in- 
homogeneous wave equation to be expressed in the 
heuristic Helmholtz integral form: 




(42) 




V 



G (xxMd^x^t {G(xxM 

(jj 





( 43 ) 



s 



8n 



- 



3G (xxM 

(JL) 



}da^ 



9n‘ 
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The first integral on the RHS of Eq. (43) contains 
the sound sources and the bulk (voliame) scattering. The 
second integral represents the surface scattering and is 
taken over all finite surfaces. 

Now, disregarding the direct path, the signal at the 
receiver is given by the surface integral alone: 



^ dn dn 



(44) 



where the subscript s denotes the scattered field. In 
order to solve Eq. (44) the following must be done; 

-Give an approximate expression for the incident 
wave . 

-Find an appropriate Greens function. 

-Make an approximation for p^(xM and 3p^/3n at 
the surface. 

As we already have assumed a simple harmonic source, 
the incident wave can be expressed as: 

(i) 



OJ 



^ ^ 2u ^ ° ^ 



R 



R 



(45) 



where 

II: power output. 

D: illumination function. 

Assuming kR>>l, i.e., that the distance of the source 
is large compared to the wavelength, the wave in the 
bounded ensonified area can be considered as a plane wave 
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characterized by its propagation vector 



k. = 
1 




thus obtaining the expression for the incident wave: 



(46) 

Further assxaming that the receiver is sufficiently far 
removed from the scattering area with the rest of the sea 
surface not contributing, then the scattering area acts as 
a small induced source in a free space and we can approxi- 
mate the propagation of the scattered waves from the en- 
sonified region in terms of the free field Greens function; 



G(X2Xi ) = 



,ikR2 



e = 



(47) 



where 




We further assume that each surface element da^ acts as a 
small reflector, and that the response of da^ to the inci- 
dent wave is that of a "local reaction," i.e., independent 
of any other part of the ensonified area A. 

Then, the Kirchhoff's approximation where it is assumed 
that p and 9p/9n vanish everywhere on the surface except at 
the ensonified area and that the values of p and 9p/9n are 
proportional to the incident wave, allow the scattered 
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target strength for the ensonified area to be approxi- 
mated as : 






= R P 



(i) 



w 



(xM 



(48) 



I = RipiiL 

3n^ A 3n^ 

where 

p is the locally reflected wave. 

R = reflection coefficient. 

_ p^c^cos 9-pc cos 9 ^ , c _ c ^ 

p^c^cos 8+pc cos 9^ sin0 sin0 ^ 

Finally, assuming a Gaussian illumination function: 




where 

X and Y are the effective dimensions of the illuminated 
area. 

The scattering integral based on the Kirchhoff's approxi- 
mation can be expressed as: 






ikBe 



ik(Ri+R2) 



2ttR 1 R 2 
,2i?(x‘y‘) 



Rf (010203) 



De 



2i (ax^+6y M 



( 50 ) 



where 

a = |- (sin 01 - sin 02 cos 83) 
S =- j sin 9 2 sin 9 3 
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Y = -^{cos 01 + cos 02) 

C(x^y') = surface displacement function. 



The above scattering integral is then applied to a randomly 
rough surface where the surface displacement function ^ is 
a random variable assumed to be represented by a Gaussian 
PDF expressed as 



Also, assume that the surface is slowly varying so that 
the signal reflects from an essentially stationary surface 
and that the succession of received scattered signals p 

n 

are assumed to form a satistically independent set of 
sample fxanctions from which sequence N the first and second 
moment of the field can be investigated. Doing so Clay and 
Tolstoy [Ref. 15] found that the mean reflected signal can 
be expressed as 




(51) 



o/2tt 



with zero mean and variance a^. 





(52) 



where 




is the signal reflected by a mirror-liJce surface, 
other factors being the same. 
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It is seen from the above that 

For o-*0 , all displacements ^ have zero proba- 
bility and the mean signal tends to Pq» Further- 
more, all elements contribute to the scattering 
coherently. 

2. For a>>k, all displacements ^ are equally probable. 
There are large phase shifts between contributions 
from different surface elements and they tend to cancel 
each other and the scattering radiation is incoherent. 

The second moment is defined as; 

T 

<s"(t) = ^ / s^(t)dt = <pp*>^ -<PP*>T (53) 



where 



sMt) = 



1 ^ 



n 



(t) - 



N<p(t)>^^ 



(54) 



In this expression the operation of squaring the signal has 
to be considered. However, each surface element has a 
different Z for a random surface. Thus, the probability 
of finding element dx' dy' with i , and element 2 with 
dx" dy" is expressed in the bivariate distribution function 
assumed to be Gaussian and of the form; 



W(ClC2)= 

where 



e[- 



(Ci^ + ?i^-2CiC2'I^) ] (55) 



2iTa^ 2(l-iJ^^)a^ 



= — <Ci(x‘y^t) ^2(x"y"t)> 
0^ 



( 56 ) 
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is the cross— correlation function characterizing 
the surface shape . 



Changing to polar coordinate leads to Eq. 6.51 of Ref. 15 
where : 



Here is the Bessel function of zero order and x is the 
transformation parameter given as 



As seen from Eq. (57) 

1. For rough surface, is large and the second 

term in the bracket, the coherent part, is 
negligible . 

2. For smooth surface, is zero and the whole 

bracket is zero. 

Since Eq. (57) cannot be integrated directly. Clay and 
Tolstoy [Ref. 15] consider it for small and large ya 
separately. 

Thus, concentrating on the high frequency limit. Clay 
and Tolstoy assumed y^a^>>! and showed that: 

-The coherent component is negligible. 

-<s^> =<pp*>, as the means tends to zero. 




cos 9 



00 




0 
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Equation (58) consists of the product of an oscillatory 
function and an exponential fxinction. Because of the Bessel 
function, the main contribution to the integral is near r=0. 
Near r=0 the phase changes slowly and the expression can be 
evaluated by the method of stationary phase. Thus, the 
expression of lii about r=0 is given as: 



Furthermore, Clay and Tolstoy show that 4^" (0) can be related 
to the characteristics of the surface as: 



Finally, Clay and Tolstoy show that the scattering signal 
can be expressed as: 




(59) 



l4^" (0) I 




(60) 



a 



>>1 





(61) 



where 



A = ensonified area 

the expected average value of Pi^, where 



Pi is the incoming pressure to the 
luminated area. 



S, ^ = scattering function 



fMe) 



e (2 y^<?^^>) 



2tt (cos 01 + cos 
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) 




'I III! ill 




f(9) — £(010203) 



1+COS0 icos 02 -sin 0 i sin0 2 cos 0 3 



COS0 1+COS02 

R = cos 0 ^ . c _ c ^ 

p^c^cos0 + pc COS0^ sin 0 sin 0 ^ 

y ~ ~2 (cos 01+ cos 02) 

K = . . 

cos 0 
k 

ct = -j (sin 01“ sin 02 cos 0 a) 




<C^^> = (3 + 5.12W) X 10 W = wind speed in m/s. 

Equation (61) is valid for: 

4y^P^ >>1 
(2 y^<C ^ ^>) >> ^ 

In summary, the reflection of high frequency signals yields 
scattered radiation which is incoherent. Furthermore, as 
pointed out by Clay and Tolstoy in [Ref. 15] although the 
radiation is primarily scattered in the specular direction, 
parts are scattered in all directions. As seen from 
Eq. (61) the scattering function is primarily dependent 

on the mean square slope of the surface and neither 

the mean square wave height nor the correlation distance 
are important. Furthermore, it is noticeable that in the 
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high frequency limit the scattering function is independent 
of the frequency since 



e 



(2y^<C^^>) 



-^1 when r^-^oo 



C. GEOMETRICAL SHADOWING 

The phenomena of shadowing of some surface areas by 
others has to be considered either when the surface irreg- 
ularities are large with respect to the incident wavelength 
or when the grazing angle is small. 

The few papers devoted to this subject are aimed mainly 
towards calculation of a "Shadowing function" based on 
the statistics of the surface. 

An explicit method, geometrical shadowing, has been 
introduced by P. Bechmann [Ref. 17] where the shadowing 
function S(9) is the probability that the point ? (Fig. 25) 
is illuminated. 

00 

S(9) = exp [-Jq{x)dx] (62) 

0 

where : 

q(x) is the probability that C is shaded by c in 
the interval (x,x+dx) given that it is not shaded 
in (0 ,x) . 

This calculation of S(9) only considers the elevation of 
the surface observation point. However, the slope also 
plays a role in that if its value exceeds cot 9 the point 
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will certainly be shaded. Thus, R. Wagner, in Ref. 18, 
incorporated both <; and using P. Bechmann's method and 

found that the conditional probability that a point on the 
surface is illuminated, given that it has height c and 
slope , can be expressed as: 

oo 

S(9i|Ci,CiM = exp[-y”q (x)dx] u(cot6 - <; M (63) 

0 

where 

u: is the unit step function. 

q(x) : is the conditional probability that c is 

shadowed in the interval (x,x+dx) given that 
it is not shadowed in (0,x) . 

The function q(x) cannot be calculated exactly. Thus, 
R. Wagner made the approximation that, for all x, the 
probability that C crosses the ray in dx is independent of 
the values of C, and at x=0. In the above, no mention 
has been made of the direction of observation. However, 
in this respect, R. Wagner [Ref. 18] points out that in 
the high frequency limit only those portions of the 
surface which are illuminated simultaneously by rays in 
the direction of incidence and observation can contribute 
to the observed scattered power. For this condition, 

R. Wagner introduced both directions as independent 
variables in the so-called bistatic shadowing probabilities 
where he defines; 
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1. S( 0 i 02 i?i/Ci') to be the conditional probability 



that the surface will not cross the incoming ray 
(Ray 1) or the outgoing ray (Ray 2) anywhere, 
given that both rays pass through an arbitrary 
point or the surface with displacement C and slope c ^ • 
2. S( 0 il 02 fCi/Ci') to be the conditional probability 
that the surface does not cross Ray 1, given that 
it does not cross Ray 2 and that both rays pass 
through the point <; having slope C ' • 

Thus, the conditional shadowing function can be expressed 



The shadowing function is then obtained by averaging 
over all possible heights and slopes 



Here w(c,C') is the bivariate PDF of the surface height 
and slope, assumed to be Gaussian 



where and ’■}>’’ are the values at x=0 of the correla- 

tion function and its second derivatives, respectively. 

For the region O<0<^ where the probability of crossing 
one ray is assumed independent of that of crossing the 
other, R. Wagner found that the bistatic shadowing function 



as : 




00 




(65) 



0 
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could be expressed as: 



{l-exp [- 2 (B 1 +B2 ) ] }x{erf Vi+erf V2} 

S(6j,02) = 

4 (B 1 +B 2 ) 



where 



B. = 
1 



exp (-V. -/if V. erfc v. 
1 1 i 

4tt V. 

1 



; i = 1,2 



V. = 






2|,,.^ti| 2Z^ 



; i = 1,2 



as we from Eq. (60) have that a^|4;"(0) | 

= <Q^ = v2 

rij_ = cot 0; i = 1,2 
and noting that 

2 -t^ 

Error functxon erf (x) = — / ® dt 

2 / -t^ 

Complementary error function erfc(x) = - J e dt 

X 

The shadowing function S(0i,02) is, in short, the fraction 
of the surface still illuminated. As seen from Eq. (60) , 
the scattered field, in the high frequency case, is pro- 
portional to the illuminated area. Hence, the shadowing 
effect of a rough surface can be introduced by multiplying 
the ensonified area A by the shadowing function S(0i,02). 
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D. ESTIMATING THE SURFACE SCATTERING EFFECT 

The following estimations are based on calculations 
in the specular direction, which is, as pointed out 
earlier, expected to give the maximum supplementary 
scattering effect. Hence, in the specular direction where 
01 = ^2 ~ ® ~ scattering fxinction Sj^^ 

reduces to the following expression: 

S . = f -- LQ L e“ ( 2 y 2 <^ i 2 >) (68) 

2 (cos 6 1 +COS 02 ) ^ 



S 



hf 



87T<i; ^ 



as 

f(9) = cose 
Y = kcos6 
a = 0 



As pointed out in the previous paragraph, the reflection 
of very high frequency signals by the sea surface yields 
scattered radiation that is incoherent under the assumption 



that 

>>1 

(2y^<C* ^>) 

Before we launch into the calculations, we will verify 
these criteria for the frequency range of interest: 60 and 

30 kHz. 
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utilizing 



cos ^8 



a = 8 . 1 X 10 



-3 



= 



aW 

46g' 



<C^^> = (3 + 5.12W) X 10 



-3 



6 = 0.74 



g = 9.81 m/s' 



and assuming a windspeed of 10 m/s (SS3) yields 
a“ = 0.30 m^ ; a = 0.55 m 



= 5.42 X 10""^ 

For 60 kHz (A = .025 m) and low grazing angles, e.g., 0 = 85° 
4y^a^ = 5.75 x 10^ >>1 



(2y^<C ^ S) = 5.2 X 10^ >> - 

and the first order 5.2>>— , R^0.5m. 

R^ 

Thus, for the 60 kHz case, the criteria are fulfilled. 
For 30 kHz (X=.05m) and 8=85°. 

= 1.44 x 10^>>! 

(2y^<C^^>) = 13>> 

R 

again the first order R>1.0 m. 

Thus, also for the 30 kHz case the criteria are fulfilled. 

Similarly, for the shadowing function in the specular 
direction 

0 < 8 i = 82 = 9 < Tr/2 and 83 = 0 
we obtain the following simplified expression 
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S{9) = 3(01,02) 



(69) 



_ [1-exp (-4B) ]erfv 
4B 



as 



Vi = V2 = V = 



"1^ (2>2)^ 



m = ri 2 = n = cot0 



Bi — B?— B— 



_ [exp (-v^ ) -/tF V erfc v] 
4 /tF V 



In suiTunary, for specular scattering at the high frequency 
limit, the expected average value of p^ , where p is the 
pressure field at the receiver, is then obtained from the 
following simplified expressions: 



<s^> = <pp*> — S 
R2 



hf 



(70) 



where 



<s^> 



= <pp*> 



AS(0) 



8irR2^<C ^ 



1-exp (-4B)erfv 

S(9) fe 



V = 



(2ZM ' 



= cot0 



B 



[exp (-v^)-/iFverfcv] 
4 (tt) V 
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~ - (3 + 5.12W) X 10 ^ ; W = windspeed in 

m/s. 



R = P^c^cosS -pc cos6^ 
P^C^COSO + P C COS0 ^ 



c__ 

sin9 



sine^ 



From Ref. 19 we use the following air/sea water inter- 
face data; 

1. For air c' = 343 m/s and = 415 Rayls. 

g 

2. For sea water c = 1500 m/s and pc = 1.54 x 10 Rayls 
The estimation of the illuminated area A for specular 

scattering where 9i = 02 = 0 and 03 = 0, is based on the 
geometry illustrated in Fig. 24b. The illuminated area 
is given by 

A = TTab (71) 

Assuming the following data to be known 
hg = depth of the source 

h = depth of the receiver 
r 

A(p = the half beam width of the directional receiver 
Both a and b of Eq. (7l) can be calculated in terms of the 
detection range R^ as follows: 



X2 




h +h 
s 



r 



( 72 ) 



Xi 



= ^ - 



X2 



= "d'I-H 



+h ' 
s r 



9 = tan ^ ( 



R, 



D 
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X 



sin9 




To a good approximation when 9 = tt/ 2, we have 



s r 



(73) 



b = R sin(A({)) 




The expression for <s^> is then introduced as a supple 
ment to the direct path to the receiver in the following 
way : 

By utilizing the relationship 




for the intensity, the scattering intensity at the receiver 
and the intensity at the ensonified area are, respectively. 



I = SL = 

S pc c 



pc pc 



Equation (70) can thus be written as 



I 



AS ( 9 ) 



s 



I 2 , , 

8ttR2 <C > 




AS ( 9 ) 



8ttR2^<c‘^> 
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lOlog I^/Ij 



lOlog = 

1 ref 



K 

o 



IL = IL(RJ = K 

5 i O 



( 75 ) 



We then have to determine IL(Ri) 



TL = SPL(l) - SPL(Rj) = 201ogRj 
SPL(Ri) = IL(Ri) = lOlogli 
SPL(l) = SL 

yielding 

IL(Ri) = SL - 201ogRi = Kj 

Then 



IL 



= K + K, 



suppl o 

is the supplement to the direct path, and 
^I-direct = SI' - ^OlogRo = K, 

, /T j. • 1 ^^direct 

Idirect/Iref ' ", 



10 



IL 



I ,/I „ = anti log -rTr = Kc 

suppl' ref ^10 ^ 

The total intensity of the receiver is thus 

, _ ^direct + ^suppl _ 

^Tot/^ref - K, + K3 



= K, 



(76) 



(77) 



(78) 



IL^^^ = lOlogKe = K? 



(79) 
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Thus, the effect of the randomly rough surface compared 
to the idealized free-field condition can be expressed as 



AIL = IL„ ^ - IL, . ,. = K„ 

Tot direct ® 



(80) 



A calculator program on a Texas Instrument 59 (later 
called TI 59) was developed to perform these calculations, 
A block diagram of the program is outlined in Fig. (26) 
and the programs steps together with a detailed description 
is given in Appendix B. The calculations are based on the 
following fixed data 
h = 2 m 



s 

h = 6 m 



r 

A(t> = 10 ° 

W = 10 m/s (SS3) 

Then, varying the detection range from R = 2000 m to 100 m 
gives the difference between IL^^^ and plotted in 

Fig. (27) . 



E. CONCLUSION AND DISCUSSION 

As seen from Fig. (27) the supplementary effect of 
the scattering from a rough surface in the high frequency 
case is negligible compared to the direct path. 

In saying so, it also should be pointed out that the 
Helmholtz-Kirchhoff ' s approach may be limited as it does 
not take into account the diffraction effects from crests 
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and throughs of the ensonified area, an effect which 
becomes increasingly important at low grazing angles, 
high frequency and when the rough surface is a super- 
position of swell and capillary waves. 
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VII. THE EFFECT OF SCATTERING AND ABSORPTION FROM THE SUB- 
SURFACE OCEAN LAYER 



A. GENERAL SCATTERING THEORY 

Generally when a region (volume) scatters sound, some 
of the energy carried by the incident wave is dispersed. 

The energy lost by the incident wave may be absorbed 
by the scatterers or it may be simply deflected from its 
original course. The amount of energy lost per second by 
the incident wave divided by the incident wave ' s intensity 
is called the total cross section of the region and is 
the sum of the absorption and scattering cross sections 




where 

Ilg = scattered power 
n = absorbed power 

The existence of gas bubbles in the subsurface ocean 
layer modifies the forward scattering in the following 
two major ways: 

1. The bubbles can resonate. When the bubbles are 

excited at a frequency near its natural frequency. 
It very efficiently absorbs and scatters the inci~ 
dent wave. At resonance, the scattering and 
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sea 



absorption cross section of a typical bubble at 
is of the order 10^ times its geometx'ical cross 
section. 

2. The bubbles change the effective compressibility of 
the water and cause the speed of sound to be a 
function of frequency, i.e., the medium is dispersive. 

We will investigate and discuss these effects by separately 
estimating: 

1. The attenuation due to the bubbles and 

2. the refraction by bubbles. 

B . ABSORPTION MODEL 

The choice of model for the subsurface ocean layer depends 
on whether the medium has a teneous or a dense distribution 
of scatterers. When the bubble density is teneous, both 
"single scattering" and "first order multiple scattering" 
approximation solutions are applicable. 

On the other extreme, when the bubble density is high, 
the so-called "diffusion" approximation can be used. 

Between these two extremes, multiple scattering effects are 
important . 

The multiple scattering theory, which in the limit also 
contains the first order approximation, will be used to 
estimate effects of attenuation due to bubbles, on the 
propagation of propeller noise from the target to the 
torpedo . 
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The geometry of the propagation model is illustrated 
in Fig. (28) where it is assumed that a plane wave is 
incident on a semi-infinite (disregarding the sea surface) 
slab of thickness x containing a number of randomly dis- 
tributed bubbles. The plane wave approximation is valid 
if the incident sound has a wavelength X much greater 
than the bubble-radius a 
ka<<l 

where 

k=cu/c=2T/X 

The receiver is located outside the slab and the beam 
pattern of the receiver is represented by the solid angle 

. 

r 

We are interested in the estimating of the total power 
received, taking into account the multiple scattering 
process in the inhomogeneous slab as well as the beam 
pattern of the receiver. 

The mathematical formulation of this problem is based 
on Twersky's theory of multiple scattering. Since the 
theory is presented in Ref. 20 only, the basic formulation, 
major assumptions, and the end results will be presented 
here . 

The total intensity is the average of the square of 
the magnitude of the total field: 
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<|4)^|^> - <|<(j;^> + ^ (81) 

where 

^ N 

'1' = 'l>j + £ u ^ 

^ s=l = 



the scalar field at the receiver location r , 

a' 

see Fig. (29) , is the sum of the incident wave (|)j^ 
and the contribution from all N scatterers. 

I I ^ is the coherent intensity based on the 
average field <|'i^£^|^> is the incoherent 

intensity based on the fluctuating field 
In Twersky's theory, the multiple scattering process is 
described by the following set of integral equations which 
Eq. (81) must satisfy: 

<U^r> = + |v <|4;2|^>p(r )d? (82) 

o & o 



where 






is an operator representing all the scattering processes 
from s to a. (See Fig. (28).) It should be noted that 
Twersky's theory includes all the multiply scattered waves 
that involve chains of successive scattering going through 
different scatterers. (See Fig. (29a).) However, the theory 
neglects the terms which include an individual 
scatterer more than once, as illustrated by Fig. (29b) . 
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Thus Twersky's theory is expected to give good results 
when back scattering is insignificant compared to the 
scattering in other directions. 

As typical for most integral equations, Twersky being 
no exception, detailed complete solutions are difficult 
to obtain. However, Twersky gave an approximate solution 
to Eq. (81) and Eq. (82) , which according to Ref. 20 have 
been found to agree reasonably well with experimental data. 
This solution is based on the following main assumptions; 

-Backscattering is assumed to be small compared to 
scattering in other directions. 

-Scattering is mostly concentrated in the forward 
direction. This is reasonable based on the assumed 
random distribution of the bubbles; i.e., no rein- 
forcement of the radiation pattern occurs except in 
the direction of the incident wave. 

-The angle 9^^ is small, i.e., 9^^ = 0°. 

This leads to the following expression for the total 
intensity at the receiver: 

<|,(;^|2> = T = exp(-po^x) exp(-pOgX) (83 

+q[l-exp(-pOgX) ] 

where 

0 ^ = absorption cross section 

0 = scattering cross section 

s 
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s 



.-n f I Mn 

4tt ' ' 



s 



where 

f = the amplitude function 

<3 = the fraction of total scattered power collected by 
the receiver as illustrated in Fig. 30 and 
P = is the bubble density, i.e., the number of 
scatterers per unit volume. 

For small values of px we see from Eq. (82) that the 
coherent part dominates: 



In this limit the multiple scattering result is equiva- 
lent to that obtained from single scattering considerations. 

For large values of px which corresponds to very dense 
or very wide slab of scatterers, the incoherent intensity 
dominates 



the receiver collects almost all the scattered power 
InT - a px 

cL 

The first case, representing the situation for teneous 
density of scatterers and/or narrow beam pattern of the 
receiver, gives a good approximation to the situation of 
interest in the thesis. It also represents the case for 



InT = -(a,+o„)px 

cL S 



(84) 



InT ~ Inq - a^px 

a 

In this latter case, it is notable that when 0^=277 



(85) 
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which no scattered power is received. This will be 
approximately true for a narrow-beam width receiver. We 
will, therefore, investigate the coherent intensity first. 

Also, the incoherent case, as represented by Eq. (86) 
will be investigated, where only losses due to absorption 
are incorporated. 



C. THE COHERENT INTENSITY CASE 

For the coherent case the intensity level after the 
incident wave has traversed a distance x is 



where 



Ix = Ip exp[- (a^+Og) px] 



Ip = incident plane wave intensity, 



(87) 



The change in intensity over the distance x is 

] 

AIL = lOlog- 



I^ I exp [-(a^+ag)px] 

i 

P P 



( 88 ) 



= lOlogexp [-(o +0 )px] 

a o 

The excess attenuation per unit distance due to bubbles 
is thus 

AIL ® 

a = - — — = - [ j 



X 



= (o^+o^) p log e 

3. S 



a = 4.34 p, in dB/m. 



(89) 



However, this only takes into account bubbles of one 
size. In a bubbly medium there is a spectr\im of radii. 
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The probability density function for finding a bubble size 

between radii a and a+da is 
- n(a) 

W(a) — ( 9 Q) 

where 

00 

/ W(a)da = 1 
o 

n(a)da is the number of bubbles per unit volume 
having radii between a and a+da. It is common to 
use da=l ym. 



p - n(a)da is the total number of bubbles per 
o 

unit volume. 

As the extinction cross section also is a function of 
the radius (See Eq. (93).), the absorption due to 
bubbles is obtained by integrating Eq. (89) over all 
possible radii 

00 

a = 4.34 J a (a)n(a)da (91) 

o 

To calculate the absorption coefficient, the extinction 
cross section must be derived from the general bubble 
dynamic relationship. This is done in detail by C. Clay 
and H. Medwin [Ref. 21] from which the follov/ing results 
are taken. 
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The scattering cross section 



a 



s 




4Tra^ 



( 92 ) 



where 



= resonance frequency = ( 

f = f = operating frequency 



2Yb3P 



A> h 



6 = damping constant “ + 6^ + 

4y 



V 



= ka + (^) (^)2 + 



Pj^wa' 



a = bubble radii 

rX(sin h X + sin X) -2 (cos h X - cos X) 



= 3(y-1) [= 



X^ (cos h X - cos X) + 3 (Y“1) (sin h X - 



3w , 

ic = a( g? 

kg = thermal conductivity of gas 
Pg = density of gas = [1 + |^-^] (1+0 . Iz) 

Pg^= density of gas at sea level 
T = surface tension 



P, = 1.013 X 10”(l+0.l2) 
A 



z = bubble depth in m. 

C = specific heat of constant pressure of gas 

pg 



zl 

sin X 
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y = shear viscosity of water 
Y = 7/5, for diatomic gas 
~ density for sea water 



b = [1 + lY.-l, sin h X - sin X j 

X COS h X - cos X 






« ■= 1 ^ - 3 ^' 



Furthermore : 



a =0+0 = 

e as 



4ira^ (6/ka) 

[(fj,/f)2 -1]2+6^ 



(93) 



and 



a =0-0 
a e s 



(94) 



A detailed computer program, as outlined in Appendix B, 
was developed for the TI 59 to handle the derivation of 
0,0, and a based on an assumed receiver depth of z = 6 m. 
For the 60 kHz case, both the extinction cross section 
a and the absorption cross section a are given in Fig. 31 
as a function of bubble radius a. 

Similarly, Fig. 32 gives o^ and for the 30 kHz case. 
Superimposed on these figures are the curves for n(a)da as 
calculated from the following: Figures 5-7 of Ref. 2 

give the resonant bubble densities in a 1 ym band as a 
function of depth and with the wind speed as parameter for 
the three discrete frequencies 12 kHz, 38 kHz and 120 kHz. 
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m 



m 



Based on these data. Fig. 33 shows the interpolated bubble 
density as a function of resonance frequency for sea state 
2 , 3 , and 6 . 

Furthermore, both A. Lovik [Ref. 2] and H. Medwin [Ref. 3] 
found that the bubble density function n(a) decreases with 
increasing bubble radii as 

-X 

naa 

where H. Medwin [Ref. 3] foimd the power law of: 

X = 4 for a<50-80 ym. 

X = 2 for a>50-80 ym. 

and A. Lovik [Ref. 2] found the power law of: 

X = 4.2 between 38 kHz and 120 kHz 
X = 2.6 between 12 kHz and 38 kHz 
which, averaged over the depth interval, corresponds to the 
bubble radii of 380 ym(12 kHz), 120 ym (38 kHz) and 49 ym 
(120 kHz) . 

As suggested by A. Lovik [Ref. 2], the discrepancy 
between the two observations is not great and may be due 
to the few measuring frequencies used in the work of 
A. Lovik . 

In summary. Fig. 33, from which we obtain the appro- 
priate resonant bubble density in a 1 ym band n(aj^) 
together with the power law naa comprise the full 
knowledge of n(a)da. 



1 
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Performing a multiplication of and n(a)da, we 
obtain Fig. 34 and Fig. 35 for the 60 kHz and the 30 kHz 
case , respectively . 

Finally, the integral 

c» 

f cr^(a)n(a)da 
0 

was evaluated using a numerical integration based on 
Simpson's discrete approximation programmed for the TI 59 
and documented in Appendix C. 

Based on the above, the following absorption coeffi- 
cient for the coherent case is obtained for f = 60 kHz 

oo 

a = 4.34 / a (a)n(a)da = 4.34 (2.016 x 10“^) 

0 

= 8.75 X 10"^ dB/m. 

a ~ 0.88 dB/m (95) 

with 

n(aj^) = 1000 

n(a) a a-*^ 
z = 6 m. 

Thus, it is seen that the attenuation due to bubbles 
is considerably greater than the "normal attenuation" 
due to chemical and viscous relaxation processes in sea 
water, which for the 60 kHz case is 
a = 0.02 dB/m. 
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Thus, the total absorption coefficient for the 60 kHz 
case is 

a ~ 0.90 dB/m. ( 95 ) 
For f = 30 kHz, the absorption coefficient due to bubbles 
in the coherent case is 



a 



= 4.34 / a (a)n(a)da = 4.34 (3.8375 x lO"^) 



with 



= 1.665 X 10 ^ dB/m. 
a ~ 0.17 dB/m 



n(aj^) = 20 



(97) 



-2 



,-l 



/ \ - 2.6 

n (a) a a 

z = 6 ra. 

For this case, the power law dependence of n(a) a a 
gives 

a = 1.655 X 10 
Thus, the difference in power law dependence makes no 
significant difference in the absorption coefficient. 

The absorption coefficient due to chemical and viscous 
relaxation processes is at 30 kHz 
a = 0.012 dB/m. 

The total absorption coefficient in the coherent 
case is 

a ~ 0.18 dB/m. (98) 
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D. THE INCOHERENT INTENSITY CASE 

For the incoherent case, where only losses due to ab- 
sorption are included, the product a and n(a)da for the 
60 kHz and 30 kHz cases are given in Fig. 36 and Fig. 37, 
respectively. Performing a numerical integration based on 
the Simpson's discrete approximation leads to the follow- 
ing results : 

For f = 60 kHz, the absorption coefficient due to 
bubbles is 

OO 1 

a = 4.34 f a (a)n(a)da = 4.34 (1.66 x 10 ) 

* 3 . 

0 

a = 0.72 dB/m. (99) 

with 

n(aj^) = 1000 

/ ^ -4 

n (a) a a 
z = 6 m. 

Adding the "normal attenuation" in sea water for 
f = kHz, yields a total absorption coefficient of 

a = 0.74 dB/m (199) 

For f = 30 kHz, the absorption coefficient due to 
bubbles is 

a = 4.34 /a (a)n(a)da = 4.34 (2.975 x lO"^) (101) 

^ 3 

0 

a = 0.13 dB/m 
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with 



n(aj^) = 20 

_2 

n(a) a a 
z = 6 m. 

Adding the "normal attenuation" in sea water for f = 30 kHz 
yields a total absorption coefficient in the incoherent case 
of 

a ~ 0.14 dB/m (102) 

E. SUMMARY AND DISCUSSION OF THE BUBBLE ATTENUATION 

As seen from the above, the attenuation due to bubbles in 
the subsurface ocean layer is important for high frequency 
and high Sea State. 

Accounting for the effect of bubbles at Sea State 3 in 
summary we found the following absorption coefficients in 
dB/m: 

f = 60 kHz f = 30 kHz 

coherent incoherent coherent incoherent 

a = 0.88 a = 0.72 a = 0.17 a = 0.13 

The main assumptions were: 

-the back scattering is small compared to that in 
other directions. 

-the scattering is mostly concentrated in the 
forward direction. 



) 
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-the angle 9^^ between the scatterer and the 

receiver is small, i.e., 9 =0°. 

as 

-the depth of the receiver is z = 6 m. 

-the sea state is 3. 

The range dependent portion of the passive sonar 
Eq. (1) TL=-201ogR-aR for both the 60 kHz and 30 kHz 
scattering results can now be compared with the reference 
data obtained in Section V. 

For f = 60 kHz, Fig. 22 gives TL=-201ogR-aR as a 
function of R with 

1. a = 0.021 dB/m, the "normal attenuation" due to 
chemical and viscous relaxation processes. 

2. a = 0.9 dB/m, the total attenuation including the 
effect of bubbles in the coherent case. 

3. a = 0.74 dB/m, the total attenuation including the 
effect of bubbles in the incoherent case. 

For the same example as in the reference model. Fig. 22 
yields the detection ranges for 60 kHz. 

R = 1200 m for a = 0.021 dB/m (103) 

R = 60 m for a = 0.90 dB/m 

R = 70 m for a = 0.74 dB/m 

Not surprisingly, this result seems to exclude the possi- 
bility or having both a searching and attack depth near 
the surface, i.e. , z = 6 m, for a torpedo system operating 
in a high frequency region, f = 60 kHz. 
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Similarly, for f = 30 kHz, Fig. 23 gives TL=-201ogR-aR 
as a function of R with 

1. a = 0.012 dB/m, the "normal attenuation" in sea 
water. 

2. ct = 0.18 dB/m, the total attenuation including the 
effect of bubbles in the coherent case. 

3. a = 0.14 dB/m, the total attenuation including the 
effect of bubbles in the incoherent case. 

For the same detection example as in the reference model. 
Fig. 23 yields the detection ranges for 30 kHz: 

R = 2400 m for a = 0.012 dB/m (104) 

R = 250 m for a = 0.18 dB/m 

R = 310 m for a = 0.14 dB/m. 

Again, the bubbles give a major decrease in the detection 
range. A detection range of R = 250 m seems marginally 
acceptable as the turn rate requirement for a pursuit 
homing trajectory may become excessive. 

The above results are summarized in Tables III and IV 
for the 60 kHz and 30 kHz, respectively. 

F. THE REFRACTION BY BUBBLES 

The presence of bubbles in the sea water affects the 
speed of sound (phase speed) primarily because of the change 
in compressibility. The derivation of this dispersive rela- 
tionship on the sound speed has been done by H. Medwin 
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DETECTION PARAMETERS AND RANGES FOR f=60 kHz 
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DETECTION PARAMETERS AND RANGES FOR f=30 kHz 
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[*R6f . 22] • H. Msdwin showed that the bubbles with resonant 
frequencies greater than the incident frequency decrease 
the sound speed, while bubbles with resonant frequency lower 
than the incident wave increase the sound speed. 

Furthermore, H. Medwin [Ref. 22] predicts the sound speed 
gradient due to bubbles as a function of depth for a frequency 
range and wind speed compatible with our domain of interest. 

He found the gradients 

g = 9c/3z = 0.26 s”^ at z = 0 m. 
g = 3c/3z = 0.016 s ^ at z = 10 m. 
g = 9c/3z = 0.005 s ^ at z = 20 ra. 

For comparison, the sound speed gradient due to pressure 
in an isothermal water is 
g = 0.017 s"^ 

This shows that the rays in the top 10 m are strongly 
influenced by bubbles. However, with respect to our surface 
scattering model (Section VI) where both the source and the 
receiver are situated very close to the surface and where 
relatively short propagation distances are encountered, 
this refraction phenomenon is assumed to have negligible 
effect . 
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VIII. THE TURN RATE LIMITATION 



As pointed out in the previous section, the presence 
of bubbles near the surface may significantly reduce the 
range at which the target can be detected. 

In this section, the turn rate necessary during 
pursuit homing at the previous estimated detection ranges 
will be studied. A computational procedure will be used 
to determine the range of angles on the bow (AOB) of the 
target at the beginning of homing which lead to miss on 
the initial attack. A trajectory where the torpedo velocity 
vector always is directed towards the instantaneous target 
position is called a pursuit homing trajectory. The 
derivation of the pursuit homing trajectory follows 
P. van Nostrand [Ref. 23] and is based on the geometry of 
Fig. 38, where 

r = approach angle, i.e., angle between ship 
velocity vector and the line of sight. 

AOB = 180 - <t> 

The equation of motion is obtained by taking the 
component along r and the normal to r, yielding 



r = Vg cos 4) - 


(105) 


r 4 ) = -V sin 4> 
s 


(106) 
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where 



r = range rate 
4) = turn rate 
Vg = target speed 

= torpedo speed. 

Dividing Eq. (105) by Eq. (106) yields: 

i Vt 1 

s 

and defining 



yields 



r 

r 



sin<j) 



- cot 4>) ^ 



(sint}))^ 

(l+COS({)) P 



Then, introducing the initial conditions: r , A 

O ' 

r^: initial detection range 



: initial approaching angle 



yields 



K = r 



(1+cos (})^) 

o"^ ! a V P -1 
(sin (p^j 



For the geometry of Fig. 38, the turning rate is 



Eq. (105) . 



(107) 

(108) 
whe re 

(109) 
given by 
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( 110 ) 



V 

0 sin $ 

Substituting Eq. (108) into Eq. (110) yields: 

P 



0 = - 



V (1+COS <P) 
s 

P-2 



( 111 ) 



K(sin 

From Eq. (108) we see that r-*-0 as 0-^0, i.e., 0 tends to 

zero as the torpedo approaches the target ship. It is 
further of interest to determine the limiting value of the 
turning rate as the torpedo approaches the target for 
various values of the parameter p. This is done by ta)cing 
the derivative of 0 with respect to 0. Thus, from Eq. (109) 
we get 



^ = - -E?(sin 0)^"^ [l+cos0]^ [2 cos - p] (112) 
a 9 i\ 

For l<p<2, we see that Eq. (112) is zero at 

cos 0 = P/2 (113) 

0 = cos ^ (P/2) 
with an associated maxima 



V (l+P/2)^[l-(P/2) 

I • I = s 

• 4>max ' K 

Furthermore, the turn rate for p<2 at impact is 
the limit of Eq. (109) yields: 

lim 0= - lim(sin 4>) ^ ^[1+cos -^0 

X\ 

0-s-O <l>-^0 



(114) 
zero as 



(115) 
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For p=2 we see from Dq. (Ill) that 



lim ^ lim(l+ cos <1)) ^ ^ (116) 

For p>2 we see from Eq. (Ill) that there is no zero 
value of d^/d(J) between (}> = 0° and $ = 180°, since all terms 
of Eq. (112) are nonzero terms for 0°<(|)<180°. Furthermore, 
the turn rate for p>2 at impact is 



I lim I 
4)^0 




(1+cos 4))^ . ^ 
P-2 ‘ 

(sin 4>) 



as (sin <J) )-*"» since (p-2)>0. 

Furthermore, as (()->-180°, we get from Eq. (112) 



(117) 



• V 

lim ^ = -j^(2+p) lim (sin 4>)^~^ (1+cos 4>)^ 

4) ->-180 

Evaluating 

lim (sin (J))^"^ (1+cos tj))^ = lim (sin 4>) ( ^gin^c)) ’ " ^ ^ 
4 >- vl 80 4)-»-180 

by applying 1' Hopitale rule to the term 



lim 



1+cos 4) 
sin (}) 



= lim 



-sin 4> 



4) -*180 



cos 4> 
4) -*180 



0 



Thus, the product sin 4 > ^ ^sin^^'^'^ approaches zero for 
any p>l, since both terms in this product approach zero 
as 43-*180. 
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Figure 39 shows a plot of computed values of |({)K/V | 
for different values of p>l. From this we can draw the 
following summarizing conclusions. 

For l<p<2 the turn rate has 

-a maximijm value at <^=cos”^ (p/2) 

-the zero value at (])=0° and <j)=180°. 

For p=2 the turn rate 
-is zero at $=180°. 

-monotonically increasing with decreasing $ 
approaching the value (4V^/K) when $=0°. 

For p 2 the turn rate 
-is zero at $=180°. 

-is monotonically increasing with decreasing $ 
approaching <» as <p=0°. 

As seen from the above for p>2, the turn rate increases 
monotonically with decreasing approaching an infinite 
turn rate to hit a point target. To avoid this singularly 
we must make some provision. If the torpedo's maximum 
turn rate is exceeded only at some very small range, a hit 
is likely. 

Figure 40 illustrates the hit criterion used. Assume 
a rudder of length L„ is situated directly behind the 
propeller which is idealizer as an acoustic point source. 

If the torpedo becomes turn rate limited at some range 
r=r', it at best can proceed along a circular path which 
lags the desired pursuit trajectory, or at worst it can 
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loose acoustic contact and go into "hold-in," maintaining 
a constant heading at the angle <{> ' • In this last case, it 
will cross the line of target ship's course in the time: 

t = r'/V^ (118) 

If this advance is less than the length of the rudder 
the torpedo will impact the rudder. Hence, the limiting 
condition is 




r- = pLj^ (119) 

For this analysis, we have arbitrarily chosen 

= 3 m (120) 

Thus, the torpedo's turn rate has not been exceeded 
when the range to the target is 

r' = pLj^ = 3 p (121) 

and a hit is assumed. 

Now we can analyze the cases p>2 and l<p<2 on the same 
footing. If we match the torpedo's maximum turn rate to 
a particular range r', it is certain that its turn rate 
has not been exceeded earlier in the run. 

This turn rate-range matching is done by substituting 
the "hit-range" defined as r' = 3 p into Eq. (106) and 
solving it for sin yielding 
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•^max 



V 

- 3 ^ sin 4 , 



sin 4 ) 




(3p) 



V 

s 



Since (p is always negative, sin 4> is always positive and 
equal to 



3pUI 



sin 4 ) = 



max 



V 



Two values of 4> satisfies this equation, 
4 »a = sin ^ ^ < ^^o 



= 180 - - 



( 122 ) 

and they are denoted 
(123) 



A unique value of K may now be found using Eq. (Ill) for 

each of the angles 4 > and 4 >„. 

A B 



K 



Vg(l+cos 4>) ^ 



These values of K are designated as K and K_, and for 

O 

each there is a corresponding value of from Eq. (109) 
with given by the initial detection range 



and ( 0.^)3 

In the case of l<p<2, the turn rate does not necessarily 
increase monotonically during the pursuit homing trajectory, 
and we must check that either 4 >a or '^’b' respectively, is not 



beyond the angle corresponding to ^ 



max 
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As seen from Eq. (113) , the turn rate reaches a maximum 
at an angle given by 

(p* = cos"^(P/2) 

with a corresponding turn rate given by 

V (l+P/2)^ [l-(P/2)^]^“^/^ 

1^1 = — = 

Thus, the value of K = K* for which the limiting turn 
rate is achieved 

_ (l+P/2) ^ [1- (P/2) 

K = 5 

l<P| 

‘ ' max 

It is important to note that, if the "A" 

trajectory is invalid since the maximum turn rate of the 
torpedo, reached at the range r = 3p at <}> is exceeded 

earlier in the trajectory. In that case, the limiting 
trajectory is the "*" trajectory. Along the same lines, 
we argue that the "B" trajectory always is a limiting 
trajectory, since the equality 4>g<<P* cannot be 

satisfied. This follows from the fact that 
< 90 ° 

and <+>g 90° as is the supplement of 

Below some critical ship speed, the torpedo will not 
be turn rate limited. This speed is obtained when the 
trajectory is normal to the ships velocity vector, <1)=90°, 
at a range of r = 3p. 
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Thus, from Eq. (106), we obtain 

sin 90° 

yielding 

(V3>nl = (3p) lil^ax 

Then introducing p = we get 

<'"s>NL = '“5) 

Furthermore, we see from Eq. (108) that for <}>=90° 

(K)^^ = 3p (126) 



Again, by probing Eq. (88) we can obtain the corres- 
ponding values of <}) , designated (<j> ),,t/ where the subscript 

O O NJ-i 

"NL" is used to indicate the "no-limit" boxindary point. 

The computational procedure is based on: 

1. Assuming a torpedo speed V = 35 kts . 

2. Using an initial detection range (beginning of 
homing) r^ = 250 m. 

3. Assuming the following maximum turn rates: 



a. 8 o/s 


b. 12 o/s 


c. 16 


o/s 


d. 24 o/s 


e. 36 o/s 


f . 48 


o/s 



4. The following range of target speeds 

0<V <25 kts 
s 

The aim of the computation is to determine whether the 
limiting and hence AOB, is governed by ({)* or by the 

turn rate at the range r = 3p. 
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The calculations are devided into two parts, and are 
performed on a TI 59 calculator. 

For a given ship speed. Part I gives the sequential 
calculations of 4>,, , <j> , K„, 0* and K* for each of the 

turn rates. The computer program is given in Appendix E. 
Then, in Part II, the probe calculation for 

~ ^‘^o^A B * performed. The program is 

given in Appendix E. 

The limiting results are given on a polar plot. Fig. 41 

with 

' 180 ° - 'Va.B,* 

As seen from Fig. 41, for a target speed of 15 kts, we 

need AOB>85° at 4 for a maximum turn rate 16 o/s in order 

o 

to have a hit at the first attack. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 



The detection performance of a passive homing torpedo 
used against shallow-draft surface ships operating in 
moderate sea states was investigated. Attention was 
focused on the effects of scattering from the randomly 
rough sea surface and scattering and absorption from the 
bubble-dominated inhomogeneous layer just below the sea 
surface. The effects of these two scattering mechanisms 
were separately estimated and their relative importance 
were compared. 

The passive sonar equation was used to predict the 
performance of the homing system, and the detection range 
considering these two scattering effects was obtained and 
compared to the detection range based on a reference 
model. An idealized propagation model was used as 
reference of comparison. This reference model was based 
on a noise source model for a cavitating propeller, the 
operational characteristics for a square-law detector, 
and a transmission model associated with a homogeneous, 
unbounded medium. 

Due to high frequency, moderate sea state and low 
grazing angles, the scattering from the randomly rough 
sea surface was found to be dominated by the direct path. 
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This result includes effects from geometrical 
shadov/ing. 

The effect of scattering and absorption from the 
bubble-dominated, inhomogeneous siibsurface layer was in- 
vestigated using multiple scattering theory. Both the 
coherent and incoherent limits were investigated by incor- 
porating the associated absorption coefficient into the 
transmission equation. 

The effect on the sound speed from the bubble content 
was found to be negligible. 

At the assumed depth setting of 6 m for the torpedo's 
search and attack phase, the scattering from the bubbles 
increased the transmission loss. This increase depended 
on the frequency and the wind speed. 

Two operating frequencies were investigated, 60 kHz 
and 30 kHz. For both cases, bubbles significantly 
decreased the detection range. 

For a torpedo system operating at the high frequencies, 
e.g., 60 kHz, the result indicates the inadvisability of 
using a searching and attack depth near the surface, i.e., 
z = 6 m. 

For an operating frequency around 30 kHz, the calculated 
detection ranges is such that the turn rate requirements 
for a pursuit homing trajectory become excessive. For a 
maximum turn rate of 16 o/s, this limitation can be avoided 
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by adapting a tactical procedure where the angle on bow 
(AOB) at the beginning of the torpedo attack is greater 
than 85°. 

At sea state 3, the results show a consistent and 
general trend towards the need for lower operating fre- 
quency in order to increase the detection range. An 
operating frequency below 30 kHz seems indicated. 

Furthermore, a search depth below the bubble-dominated 
subsurface layer, i.e., z>15-20 m would result in an 
increased detection range. 

To reduce the operational limitations induced by the 
scattering and absorption effects, a high maximum turn 
rate together with a variable speed capability, where 
p-2 would be beneficial. 

The result of this analysis has clearly demonstrated 
the importance of environmental factors on the torpedo 
capability, and is useful in giving insight into the be- 
havior of a homing torpedo during its critical attack 
phase . 

A valuable follow-on of this study would be an investi- 
gation of the effects of the bubble-dominated subsurface 
layer on target validation and pitch plane steering when 
the torpedo search depth is 50-60 m. 
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APPENDIX A 



DETAILED OCEANOGRAPHIC BACKGROUND MATERIAL 

A. GEOGRAPHIC DESCRIPTION 

The Norwegian coastal waters constitute the eastern 
boundary of the Norwegian Sea. Although some general 
aspects related to the Norwegian Sea will be covered, 
this analysis will be concentrated on the Norwegian 
coastal waters above 68°N. 

B. WIND 

The northern region of the Norwegian Sea is affected 
by the Polar Easterlies and the southern region by the 
prevailing Westerlies. There are two dominant air masses 
which are relatively permanent; 

-The Greenland high. 

-The Iceland low. 

These pressure systems produce storms which are carried 
across the Norwegian Sea in a belt from Iceland towards 
the Norwegian Coast causing steady precipitation and wind 
most of the year. The steep Norwegian Coast has a con- 
siderable influence on the winds and consequently also on 
the waves in the coastal waters. The main general modi- 
fications are that the streamlines tend to run parallel 
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to the coast and that wind and sea increases with distcince 
from the sheltered coast into open ocean. 

Strong local variation may occur. The most important 
of these are the marked local increase in wind speed in 
areas where the coast sharply changes direction. One such 
"corner effect," caused by the confluence of the streamlines, 
occurs near "Nordkapp" (North Cape) . Also of importance 
are monsoonal effects due to the different heat capacities 
between the ocean and the continent. Drainage of cold air 
from the inland valleys in the wintertime causes a marked 
increase in the wind speed in several areas along the 
coast. Most of these coastal effects are significantly 
dissipated at distances of approximately 50 nmi. from the 
coast. The fact that the wind tends to blow along the 
coast is clearly demonstrated in Figs. 12 and 13 which 
include data from weather stations from "Hillesoy" to 
"Ona" and "Myken" to "Furuholmen" respectively. The high 
frequency of offshore winds is caused by the drainage of 
cold air from the inland valleys during winter time. 

That this phenomenon is closely connected to the coast 
is illustrated by the fact that it is missing at weather 
station "Skomvaer" situated approximately 50 nmi. off the 
main coast. A frequency distribution of observed wind 
speeds along the coast, obtained from Ref. 1, is presented 
in Table V. A summary of this table follows: 
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TABLE V 



FREQUENCY DISTRIBUTION OF WIND SPEED IN PERCENT PER YEAR 
AT WEATHER STATIONS ALONG THE NORWEGIAN COAST 



m/sec 


0-7 


8-13 


14-20 


^21 


Beaufort 


0-4 


5+6 


7+8 


=9 


Ferder 


67 


28 


5 


0.1 


Lyngor 


79 


19 


2 


0.2 


Lista 


67 


28 


5 


0.2 


Utsira 


73 


22 


5 


0.3 


Hellisoy 


77 


19 


4 


0.2 


Krakenes 


58 


28 


11 


2.7 Jan-Dec 










^ ^ 1949-1975 


Ona 


70 


22 


7 


0.7 


Sula 


61 


30 


8 


1.1 


Nordoyan 


49 


35 


14 


1.9 


Myken 


64 


26 


9 


1.2 


Skomvaer 


58 


31 


10 


0.9 


Andenes 


79 


18 


3 


1 — 1 

• 

o 


Torsvar 


72 


23 


5 


0.4 


Fruholmen 


54 


32 


12 


1.8 


Vardo 


75 


22 


3 


0.1 


Biornoya 


63 


31 


6 


0.4 1956-1975 


Polarfront 


46 


30 


14 


1.4 
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The highest winds are reported from the areas between 
62 N and 68*^N. The frequency of high winds in this area 
are significantly greater than those found at stations to 
the south of 62°N and also to the north of 68°N. The wind 
condition around North Cape, are very severe. This area 
can be compared with the other coastal area of high wind 
speed such as "Stadt" (represented by the weather station, 
"Krakenes"). In these areas the frequency of storms is 
greater than at "Polarfront” situated in the open Norwegian 
Sea. 



C . WAVES 

Frequency distribution of significant wave heights are 
represented in Table VI. 

The station north of 68°N is characterized by compara- 
tively small frequency of high waves. Even at "Furuholmen" 
where wind conditions are very severe, the frequency dis- 
tribution of significant wave heights is similar to more 
sheltered areas like "Utsira." 

The seasonal variation for the area of interest around 
70°N is given in Figs. 3 and 4. The average monthly dis- 
tribution of significant wave heights for a typical 
station like "Andenes" is given in Fig. 3 and Fig. 4 
gives the yearly distribution for this station. 
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TABLE VI 



FREQUENCY DISTRIBUTION OF SIGNIFICANT WAVE 
HEIGHT IN PERCENT AT THE WEATHER STATION "ANDENES" 



STATION: Andenes YEARS: 1949-1972 



SS 


0+1+2 


3 


4 


5 


6 


7 


8 


9 




^s 


0 


0.5 - 


1.3 


- 2.5 - 


4.0 


6.0 


- 9.0- 


13.9 


N 




0.4 


1.2 


2.4 


3.9 


5. 9 


8.9 


13.9 






JAN 


15.3 


26.9 


33.3 


15.5 


6.6 


2.3 


0.18 


0.00 


2542 


FEB 


17. 7 


28.7 


31.3 


14.2 


5.9 


2.0 


0.14 


0.00 


2295 


MAR 


18.1 


28.9 


30.8 


14.3 


5.8 


1.9 


0.17 


0.00 


2529 


APR 


24.5 


33.1 


27.4 


10.8 


3.4 


0.8 


0.04 




2293 


lAAY 


30.9 


35.5 


23.9 


7.7 


1.7 


0.3 


0.01 




1720 


JUN 


37.4 


36.9 


20.2 


4.8 


0.6 


0.1 


0.00 




1565 


JUL 


42.6 


36.2 


17.1 


3.7 


0.4 


0.0 


0.00 




1618 


AUG 


40.6 


36.4 


18.0 


4.3 


0.5 


0.1 






1621 


SEP 


29.7 


34.3 


24.4 


8.9 


2.1 


0.5 


0.01 




2301 


OCT 


21.8 


31.3 


29.4 


12.4 


4.1 


1.0 


0.01 




2374 


NOV 


21.7 


31. 5 


29.9 


11.7 


4.1 


1.1 


0.06 




2303 


DEC 


18.5 


30.1 


32.2 


13.2 


4.7 


1.2 


0.09 




2379 


YEAR 


26.6 


32.5 


26.5 


10.1 


3.3 


0.9 


0.06 


0.00 





N = Number of Observations 
SS = Class Interval, State of Sea 
Hg = Significant Wave Height 
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D. AMBIENT NOISE 



Due to the fact that the Norwegian Sea is physically 
separated from the Atlantic by the Faeroy-Shetland-Iceland 
ridge, and from the Greenland Sea by the Jan Mayan ridge, 
long distance shipping noise is transferred into 
the area. This, combined with relatively low shipping 
traffic in the central and northern parts of the Norwegian 
Sea, produces a relatively low ambient noise level for the 
frequency band 100<f<1000 Hz, especially when noise from 
marine life is not included. Further, in the absence of 
nearby shipping and marine life, the ambient noise level 
in the frequency band Kf<50 kHz is, according to Ref. 4, 
dominated by the wind. 

A typical area for the central part of the Norwegian 
Sea can be represented by the weather station "Polarfront" 
at 66°N, 2°E. Reference 1 shows that there is approximately 
a 15% chance of finding wind forces of Beaufort>6. These 
effects predict a moderate ambient noise level in the 
frequency range 100 Hz-50 kHz for the central and northern 
regions of the Norwegian Sea. 

The shallow coastal Norwegian waters are, according to 
Ref. 4, typically 5-10 dB noisier than the corresponding 
deep waters. However, great variability caused by local 
traffic, fishing fleet activity, marine life, and local 
wind conditions makes ambient noise level prediction 
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difficult in these areas. This means that accurate ambient 
noise level determinations have to be made on the spot, as 
it is both site and time dependent. 

A noticeable influence on the ambient noise level is 
rain, which is a year around feature along the Norwegian 
coast. As seen from Fig. 8, taken from Ref. 4, rain has 
a tendency to produce a constant high ambient noise level 
over a large frequency range, thus dominating other 
effects. Furthermore, for the upper frequency of interest, 
i.e., around 60 kHz the lower bound for the ambient 
noise is determined by the thermal agitation, see Fig. 8. 

In determining the figure of merit (FOM) for a 
passive sonar system, the noise level will be the larger 
of either the self noise or the ambient noise. For a 
torpedo the self noise will typically be dominant. 

E. SOUND SPEED PROFILES 

Again, concentrate on data relevant to Norwegian 
coastal waters. According to Ref. 6, which covers the 
southern part of the Norwegian coast, low sound speeds 
are common because of the influence of water from the 
Baltic Sea combined with fresh water drainage from the 
fjords. Furthermore, great variability, both seasonal 
and within seasons, is encountered. Figure 9 obtained 
from Ref. 5 gives a picture of the sound speed profiles 
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for the northern Norwegian coast. Again, large variations 
are common. Noticeable in both sets of data is a typical 
seasonal pattern of strong cooling of the surface layer 
during winter and a similarly strong heating during 
summer. Furthermore, note that the minimum and maximum 
are relatively shallow, i.e., less than 50 m. 

Also characteristic is the influence of the cold and 
fresh melt waters drained out through the fjord-arms 
during spring and summer. 

To illustrate the sonar problems associated with these 
sound speed profiles, ray paths for the extremes of Fig. 9 
are shown in Figs. 10 and 11, where the source is 3 m below 
the sea surface. 
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APPENDIX B 



SURFACE SCATTERING TI 59 PROGRAM 
A. INTRODUCTION 

This program gives specular scattered power at the 
receiver versus incoming power at the randomly rough 
surface in the high frequency limit according to Eq. (70) . 

Shadowing of surface areas by other parts of the 
boundary are taken into accoiint by the bistatic shadowing 
function S ( 9) . 

Furthermore, the program gives the effect of the 
randomly rough surface compared to the idealized free- 
field condition as expressed in Eq. (80) . The results 
of these calculations is given in Fig. 27. 

3 . PROGRAM STEPS 

A block-diagram of the computer program is given in 
Fig. 26. 

The program uses the partitioning ratio of program 
to data space according to code 4 OP17. The users in- 
structions are as follows : 



no 



Procedure 


Enter 


Press 


Display 


Enter data 


Detection range 


2nd A 


1 — I 

P 


Enter data 


Source depth 


R/S 


h 1 
s 


Enter data 


Receiver depth 


R/S 


h 1 
r 


Enter data 


Beam width 


R/S 




Enter data 


Wind speed 


R/S 




Calculate 9, 


«... 


2nd B 


A 


/ Ra / A 






(0,Ri,R2,A) 


Calculate 


— 


2nd C 




Calculate v 


— 


2nd D 


vl 


Calculate erfc v 


— 


A 


erfc v^ 


Calculate S(9) 


-- 


B 


S(0)^ 


<s^> 






<s^>^ 


Calculate <PiPi*> 


— 


C 


A 

•— * 

V 


Calculate AIL 


— 


D 


AIL^ 



^These values are printed automatically if the calculator 
is connected to the PC~100A Print Cradle. 



For the error function complement we have 



erfc(v) = — y e ^ da 

/tt V 

in contrast to the normal distribution 

a(u) = 

/2n a 

However, there is a linear relationship between the 



two functions. 



Ill 



The numerical equation used to calculate the erfc(v) 
is a modified program from Texas Instruments [Ref. 24]. 
erfc(v) = Z (v) [bia+b2a^+b3a^+b4a‘*+b5a5] 

where 



p = .231649 

b = .451673691 
1 

b2 = -.504257336 
b3 = 2.51939026 
b^ = -2.563346623 
bs = 1.881292139 



The program steps are listed below, giving location 



(LOG) , code (COD) , key symbol (KEY) , and comments. 



000 


76 


LBL 


023 


17 


B' 


046 


05 


05 


001 


16 


A' 


024 


53 


( 


047 


55 


■r 


002 


42 


ST 


025 


53 


( 


048 


43 


RCL 


003 


00 


00 


026 


43 


RCL 


049 


02 


02 


004 


99 


PRT 


027 


00 


00 


050 


54 


) 


005 


91 


R/S 


028 


65 


X 


051 


22 


INV 


006 


42 


STO 


029 


43 


RCL 


052 


30 


TAN 


007 


01 


01 


030 


02 


02 


053 


54 


) 


008 


99 


PRT 


031 


54 


) 


054 


42 


STO 


009 


91 


R/S 


032 


55 




055 


06 


06 


010 


42 


STO 


033 


53 


( 


056 


99 


PRT 


Oil 


02 


02 


034 


43 


RCL 


057 


53 


( 


012 


99 


PRT 


035 


01 


01 


058 


53 


( 


013 


91 


R/S 


036 


85 


+ 


059 


43 


RCL 


014 


42 


STO 


037 


43 


RCL 


060 


05 


05 


015 


03 


03 


038 


02 


02 


061 


55 


V 


016 


99 


PRT 


039 


54 


) 


062 


53 


( 


017 


91 


R/S 


040 


54 


) 


063 


43 


RCL 


018 


42 


STO 


041 


42 


STO 


064 


06 


06 


019 


04 


04 


042 


05 


05 


065 


38 


SIN 


020 


99 


PRT 


043 


53 


( 


066 


54 


) 


021 


91 


R/S 


044 


53 


( 


067 


54 


) 


022 


76 


LBL 


045 


43 


RCL 


068 


42 


STO 



112 



069 

070 

071 

072 

073 

074 

075 

076 

077 

078 

079 

080 

081 

082 

083 

084 

085 

086 

087 

088 

089 

090 

091 

092 

093 

094 

095 

096 

097 

098 

099 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 



07 


07 


120 


43 


RCL 


171 


69 


OP 


99 


PRT 


121 


03 


03 


172 


06 


06 


53 


( 


122 


38 


SIN 


173 


91 


R/S 


53 


( 


123 


65 


X 


174 


76 


LBL 


43 


RCL 


124 


43 


RCL 


175 


19 


D' 


00 


00 


125 


07 


07 


176 


53 


( 


75 


- 


126 


54 


) 


177 


43 


RCL 


43 


RCL 


127 


54 


) 


178 


06 


06 


05 


05 


128 


42 


STO 


179 


30 


TAN 


54 


) 


129 


10 


10 


180 


35 


1/x 


55 


-f 


130 


69 


OP 


181 


54 


) 


53 


( 


131 


06 


06 


182 


42 


STO 


43 


RCL 


132 


91 


R/S 


183 


12 


12 


06 


06 


133 


76 


LBL 


184 


53 


( 


38 


SIN 


134 


18 


C ' 


185 


43 


RCL 


54 


) 


135 


25 


CLR 


186 


12 


12 


54 


) 


136 


69 


OP 


187 


55 


•=• 


42 


STO 


137 


00 


00 


188 


53 


( 


08 


08 


138 


00 


00 


189 


53 


( 


99 


PRT 


139 


06 


6 


190 


02 


2 


53 


( 


140 


02 


2 


191 


65 


X 


43 


RCL 


141 


07 


7 


192 


43 


RCL 


07 


07 


142 


03 


3 


193 


11 


11 


85 


+ 


143 


02 


2 


194 


54 




43 


RCL 


144 


03 


3 


195 


34 




08 


08 


145 


03 


3 


196 


54 


) 


54 


) 


146 


69 


OP 


197 


54 


) 


42 


STO 


147 


04 


04 


198 


42 


STO 


09 


09 


148 


53 


( 


199 


13 


13 


25 


CLR 


149 


53 


( 


200 


99 


PRT 


69 


OP 


150 


53 


( 


201 


91 


R/S 


00 


00 


151 


43 


RCL 


202 


76 


LBL 


01 


1 


152 


04 


04 


203 


11 


A 


03 


3 


153 


65 


X 


204 


53 


( 


03 


3 


154 


05 


5 


205 


53 


( 


05 


5 


155 


93 


. 


206 


53 


( 


01 


1 


156 


01 


1 


207 


53 


( 


07 


7 


157 


02 


2 


208 


43 


RCL 


01 


1 


158 


54 


) 


209 


13 


13 


03 


3 


159 


85 


+ 


210 


33 


x^ 


69 


OP 


160 


03 


3 


211 


54 


( 


04 


04 


161 


54 


) 


212 


55 


V 


53 


( 


162 


65 


X 


213 


02 


2 


53 


( 


163 


93 


• 


214 


54 


) 


89 


TT 


164 


00 


0 


215 


94 


+/- 


65 


X 


165 


00 


0 


216 


22 


INV 


43 


RCL 


166 


01 


1 


217 


23 


LNX 


08 


08 


167 


54 


) 


218 


54 


) 


54 


) 


168 


50 


ixl 


219 


65 


X 


65 


X 


169 


42 


STO 


220 


53 


( 


53 


( 


170 


11 


11 


221 


53 


( 



113 



222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 



89 


II 


273 


03 


3 


324 


43 


RCL 


34 




274 


08 


8 


325 


15 


15 


65 


X 


275 


01 


1 


326 


45 


YX 


02 


2 


276 


05 


5 


327 


04 


4 


54 


) 


277 


03 


3 


328 


85 


+ 


35 


1/x 


278 


75 


- 


329 


01 


1 


54 


) 


279 


93 


• 


330 


93 


• 


54 


) 


280 


03 


3 


331 


03 


3 


42 


STO 


281 


05 


5 


332 


03 


3 


14 


14 


282 


06 


6 


333 


00 


0 


53 


( 


283 


05 


5 


334 


02 


2 


53 


( 


284 


06 


6 


335 


07 


7 


53 


( 


285 


03 


3 


336 


04 


4 


53 


( 


286 


07 


7 


337 


04 


4 


43 


RCL 


287 


08 


8 


338 


02 


2 


13 


13 


288 


02 


2 


339 


09 


9 


65 


X 


289 


65 


X 


340 


00 


0 


93 


• 


290 


43 


RCL 


341 


65 


X 


02 


2 


291 


15 


15 


342 


43 


RCL 


03 


3 


292 


33 


x^ 


343 


15 


15 


01 


1 


293 


85 


+ 


344 


45 


yX 


06 


6 


294 


01 


1 


345 


05 


5 


04 


4 


295 


93 


• 


346 


54 


) 


01 


1 


296 


97 


7 


347 


54 


) 


09 


9 


297 


08 


8 


348 


42 


STO 


86 


+ 


298 


01 


1 


349 


16 


16 


01 


1 


299 


04 


4 


350 


99 


PRT 


54 


) 


300 


07 


7 


351 


91 


R/S 


35 


1/x 


301 


07 


7 


352 


76 


LBL 


54 


) 


302 


09 


9 


353 


12 


B 


42 


STO 


303 


03 


3 


354 


53 


( 


15 


15 


304 


07 


7 


355 


53 


( 


53 


( 


305 


65 


X 


356 


53 


( 


53 


( 


306 


43 


RCL 


357 


43 


RCL 


43 


RCL 


307 


15 


15 


358 


13 


13 


14 


14 


308 


45 


YX 


359 


33 


x^ 


65 


x 


309 


03 


3 


360 


94 


+/- 


01 


1 


310 


75 


- 


361 


22 


INV 


93 




311 


01 


1 


362 


23 


LNX 


08 


8 


312 


93 


• 


363 


54 


) 


01 


1 


313 


08 


8 


364 


75 


- 


54 


) 


314 


02 


2 


365 


53 


( 


65 


X 


315 


01 


1 


366 


89 


IT 


53 


( 


316 


02 


2 


367 


34 


/x 


43 


RCL 


317 


05 


5 


368 


65 


X 


15 


15 


318 


05 


5 


369 


43 


RCL 


65 


X 


319 


09 


9 


370 


13 


13 


93 




320 


07 


7 


371 


65 


X 


03 


3 


321 


08 


8 


372 


43 


RCL 


01 


1 


322 


00 


0 


373 


16 


16 


09 


9 


323 


65 


X 


374 


54 


) 



114 



375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 



54 


) 


426 


99 


PRT 


477 


14 


D 


55 


T 


427 


91 


R/S 


478 


25 


CLR 


53 


( 


428 


76 


LBL 


479 


69 


OP 


89 


TT 


429 


13 


C 


480 


00 


00 


34 




430 


25 


CLR 


481 


01 


1 


65 


X 


431 


69 


OP 


482 


06 


6 


04 


4 


432 


00 


00 


483 


01 


1 


65 


X 


433 


03 


3 


484 


04 


4 


43 


RCL 


434 


06 


6 


485 


69 


OP 


13 


13 


435 


01 


1 


486 


04 


04 


54 


) 


436 


05 


5 


487 


53 


( 


54 


) 


437 


01 


1 


488 


53 


( 


42 


STO 


438 


03 


3 


489 


43 


RCL 


17 


17 


439 


03 


3 


490 


19 


19 


53 


( 


440 


07 


7 


491 


28 


LOG 


53 


( 


441 


69 


OP 


492 


54 


) 


53 


( 


442 


04 


04 


493 


65 


X 


01 


1 


443 


53 


( 


494 


01 


1 


75 


— 


444 


53 


{ 


495 


00 


0 


53 


( 


445 


43 


RCL 


496 


54 


) 


53 


{ 


446 


10 


10 


497 


42 


STO 


53 


( 


447 


65 


X 


498 


20 


20 


43 


RCL 


448 


43 


RCL 


499 


69 


OP 


17 


17 


449 


18 


18 


500 


06 


06 


65 


X 


450 


54 


) 


501 


91 


R/S 


04 


4 


451 


55 


T 








54 


) 


452 


53 


( 








94 


+/- 


453 


53 


( 








54 


) 


454 


53 


( 








22 


IW 


455 


53 


( 








23 


LNX 


456 


43 


RCL 








54 


) 


457 


07 


07 








54 


) 


458 


33 


x^ 








65 


X 


459 


54 


) 








53 


( 


460 


65 


X 








01 


1 


461 


08 


8 








75 


— 


462 


54 


) 








43 


RCL 


463 


65 


X 








16 


16 


464 


89 


IT 








54 


) 


465 


54 


) 








54 


) 


466 


65 


X 








55 




467 


43 


RCL 








53 


( 


468 


11 


11 








43 


RCL 


469 


54 


) 








17 


17 


470 


54 


) 








65 


X 


471 


42 


STO 








04 


4 


472 


19 


19 








54 


) 


473 


69 


OP 








54 


) 


474 


06 


06 








42 


STO 


475 


91 


R/S 








18 


18 


476 


76 


LBL 
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APPENDIX C 



BUBBLE DYNAMICS TI 59 PROGRAM 

A. INTRODUCTION 

This program gives the resonant frequency and the 
damping constant 6 for bubbles according to Eq. (92) . 
Furthermore, the program gives scattering cross section 
0 g, extinction cross section and absorption cross 
section as a function of bubble radius a, incoming fre- 
quency f, and depth z. The results of the calculations 
are given in Figs. 31 and 32. 

B . PROGRAM STEPS 



The user's instructions are as follows: 

Procedure Enter Press Display 



Enter data 


Bubble Radius 


2nd A 


a^ 


Enter data 


Incoming Frequency 


R/S 


f^ 


Enter data 


Depth 


R/S 




Calculate X 


— 


2nd B 


X^ 


Calculate 
coshX and 
sinhX 




2nd C 


coshX^ 

sinhX^ 


Calculate d/b 


— 


2nd D 


d/b^ 


Calculate b 


— 


2nd E 


b^ 


Calculate 6 


— 


A 




Calculate f^ 


— 


B 


f 1 
r 


Calculate 5 


— 


C 




Calculate a 


— 


D 


0 ^ 
e 


Calculate o 
and a 


— 


E 


0 1 
s 
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^These values are printed automatically if the calculator 
is connected to PC-IOOA Print Cradle. 



The program steps are listed below giving location (LOG) 



code (COD) , key (KEY) and comments. 



000 


76 


LBL 


043 


55 


T 


086 


01 


01 


001 


16 


A ' 


044 


43 


RCL 


087 


99 


PRT 


002 


42 


STO 


045 


11 


11 


088 


91 


R/S 


003 


11 


11 


046 


54 


) 


089 


76 


LBL 


004 


99 


PRT 


047 


85 


+ 


090 


18 


C ’ 


005 


91 


R/S 


048 


93 


• 


091 


53 


( 


006 


42 


STO 


049 


00 


0 


092 


53 


( 


007 


12 


12 


050 


00 


0 


093 


43 


RCL 


008 


99 


PRT 


051 


01 


1 


094 


01 


01 


009 


91 


R/S 


052 


02 


2 


095 


22 


INV 


010 


42 


STO 


053 


09 


9 


096 


23 


LNX 


Oil 


13 


13 


054 


65 


X 


097 


85 


+ 


012 


99 


PRT 


055 


53 


( 


098 


43 


RCL 


013 


91 


R/S 


056 


01 


1 


099 


01 


01 


014 


76 


LBL 


057 


85 


+ 


100 


94 


X/- 


015 


17 


B ' 


058 


93 


m 


101 


22 


INV 


016 


53 


( 


059 


01 


1 


102 


23 


LNX 


017 


53 


( 


060 


65 


X 


103 


54 


) 


018 


53 


( 


061 


43 


RCL 


104 


55 




019 


53 


( 


062 


13 


13 


105 


02 


2 


020 


53 


( 


063 


54 


) 


106 


54 


) 


021 


53 


( 


064 


54 


) 


107 


42 


STO 


022 


04 


4 


065 


54 


) 


108 


02 


02 


023 


65 


X 


066 


65 


X 


109 


99 


PRT 


024 


89 


TT 


067 


93 




110 


53 


( 


025 


54 


) 


068 


02 


2 


111 


53 


( 


026 


65 


X 


069 


04 


4 


112 


43 


RCL 


027 


43 


RCL 


070 


54 


) 


113 


01 


01 


028 


12 


12 


071 


55 


T 


114 


22 


INV 


029 


54 


) 


072 


93 


• 


115 


23 


LNX 


030 


65 


X 


073 


00 


0 


116 


75 


— 


031 


53 


( 


074 


00 


0 


117 


43 


RCL 


032 


53 


( 


075 


00 


0 


118 


01 


01 


033 


93 




076 


00 


0 


119 


94 


+/- 


034 


00 


0 


077 


05 


5 


120 


22 


INV 


035 


00 


0 


078 


06 


6 


121 


23 


LNX 


036 


00 


0 


079 


54 


) 


122 


54 


) 


037 


00 


0 


080 


34 


/x 


123 


55 


T 


038 


00 


0 


081 


65 


X 


124 


02 


2 


039 


00 


0 


082 


43 


RCL 


125 


54 


) 


040 


01 


0 


083 


11 


11 


126 


42 


STO 


041 


09 


9 


084 


54 


) 


127 


03 


03 


042 


01 


1 


085 


42 


STO 


128 


99 


PRT 



117 



129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 



91 


R/S 


180 


53 


( 


231 


54 


) 


76 


LBL 


181 


01 


1 


232 


54 


) 


17 


D' 


182 


93 


. 


233 


55 




53 


( 


183 


02 


2 


234 


53 


( 


53 


( 


184 


65 


X 


235 


43 


RCL 


53 


( 


185 


43 


RCL 


236 


02 


02 


53 


( 


186 


01 


01 


237 


75 


— 


43 


RCL 


187 


54 


) 


238 


43 


RCL 


01 


01 


188 


65 


X 


239 


01 


01 


65 


X 


189 


53 


( 


240 


39 


COS 


53 


{ 


190 


43 


RCL 


241 


54 


} 


43 


RCL 


191 


03 


03 


242 


54 


) 


03 


03 


192 


75 


- 


243 


85 


+ 


85 


+ 


193 


43 


RCL 


244 


01 


1 


43 


RCL 


194 


01 01 


245 


54 


) 


01 


01 


195 


38 


SIN 


246 


35 


1/x 


38 


SIN 


196 


54 


} 


247 


65 


X 


54 


) 


197 


54 


) 


248 


53 


( 


54 


} 


198 


54 


} 


249 


54 


( 


75 


- 


199 


54 


) 


250 


43 


RCL 


02 


2 


200 


65 


X 


251 


04 


04 


65 


X 


201 


01 


1 


252 


33 


x2 


53 


( 


202 


93 


• 


253 


54 


) 


43 


RCL 


203 


02 


2 


254 


85 


+ 


02 


02 


204 


54 


) 


255 


01 


1 


75 


- 


205 


42 


STO 


256 


54 


) 


43 


RCL 


206 


04 


04 


257 


35 


1/x 


01 


01 


207 


99 


PRT 


258 


54 


) 


39 


COS 


208 


91 


R/S 


259 


42 


STO 


54 


) 


209 


76 


LBL 


260 


05 


05 


54 


} 


210 


10 


E 


261 


99 


PRT 


55 


T 


211 


53 


( 


262 


91 


R/S 


53 


( 


212 


53 


( 


263 


76 


LBL 


53 


( 


213 


53 


( 


264 


11 


A 


53 


( 


214 


53 


( 


265 


53 


( 


43 


RCL 


215 


53 


( 


266 


53 


( 


01 


01 


216 


03 


3 


267 


53 


( 


33 


X2 


217 


93 


. 


268 


01 


1 


54 


) 


218 


02 


2 


269 


05 


5 


65 


X 


219 


55 


-t. 


270 


00 


0 


53 


( 


220 


43 


RCL 


271 


55 


T 


43 


RCL 


221 


01 


01 


272 


53 


( 


02 


02 


222 


54 


) 


273 


01 


1 


75 


- 


223 


65 


X 


274 


00 


0 


43 


RCL 


224 


53 


( 


275 


01 


1 


01 


01 


225 


43 


RCL 


276 


03 


3 


39 


COS 


226 


03 


03 


277 


00 


0 


54 


) 


227 


75 


- 


278 


00 


0 


54 


) 


228 


43 


RCL 


279 


00 


0 


85 


+ 


229 


01 


01 


280 


65 


X 


53 


( 


230 


38 


SIN 


281 


43 


RCL 
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1l 

















282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 



11 


11 


333 


06 


06 


384 


53 


( 


65 


X 


334 


65 


X 


385 


93 




53 


( 


335 


01 


1 


386 


00 


0 


01 


1 


336 


00 


0 


387 


00 


0 


85 


+ 


337 


01 


1 


388 


00 


0 


93 


• 


338 


03 


3 


389 


00 


0 


01 


1 


339 


00 


0 


390 


04 


4 


65 


X 


340 


00 


0 


391 


01 


1 


43 


RCL 


341 


00 


0 


392 


08 


8 


13 


13 


342 


54 


) 


393 


08 


8 


54 


) 


343 


65 


X 


394 


07 


7 


54 


) 


344 


53 


( 


395 


09 


9 


54 


) 


345 


01 


1 


396 


65 


X 


65 


X 


346 


85 


+ 


397 


43 


RCL 


53 


) 


347 


93 


• 


398 


12 


12 


01 


1 


348 


01 


1 


399 


65 


X 


75 


— 


349 


65 


X 


400 


43 


RCL 


53 


( 


350 


43 


RCL 


401 


11 


11 


53 


( 


351 


13 


13 


402 


54 


) 


04 


4 


352 


52 


) 


403 


85 


+ 


93 


• 


353 


54 


) 


404 


53 


( 


02 


2 


354 


55 




405 


43 


RCL 


65 


X 


355 


01 


1 


406 


04 


04 


43 


RCL 


356 


93 


• 


407 


65 


X 


05 


05 


357 


00 


0 


408 


53 


( 


54 


) 


358 


02 


2 


409 


53 


( 


35 


1/x 


359 


06 


6 


410 


43 


RCL 


54 


) 


360 


54 


) 


411 


07 


07 


54 


) 


361 


50 


1x1 


412 


55 




54 


) 


362 


34 




413 


43 


RCL 


85 


+ 


363 


65 


X 


414 


12 


12 


01 


1 


364 


53 


( 


415 


54 


) 


54 


) 


365 


53 


( 


416 


33 


X2 


42 


STO 


366 


89 


TT 


417 


54 


) 


06 


06 


367 


65 


X 


418 


54 


) 


99 


PRT 


368 


02 


2 


419 


54 


) 


91 


R/S 


369 


65 


X 


420 


85 


+ 


76 


LBL 


370 


43 


RCL 


421 


53 


( 


12 


B 


371 


11 


11 


422 


93 


= 


53 


( 


372 


54 


) 


423 


00 


0 


53 


( 


373 


35 


1/x 


424 


00 


0 


53 


( 


374 


54 


) 


425 


06 


6 


53 


( 


375 


54 


) 


426 


02 


2 


04 


4 


376 


42 


STO 


427 


00 


0 


93 


. 


377 


07 


07 


428 


05 


5 


02 


2 


378 


99 


PRT 


429 


55 


T 


65 


X 


379 


91 


R/S 


430 


53 


( 


43 


RCL 


380 


76 


LBL 


431 


43 


RCL 


05 


05 


381 


13 


C 


432 


12 


12 


65 


X 


382 


53 


( 


433 


65 


X 


43 


RCL 


383 


53 


( 


434 


53 


( 
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435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 



43 


RCL 


486 


43 


RCL 


537 


33 


x2 


11 


11 


487 


07 


07 


538 


54 


) 


33 




488 


55 


■r 


539 


75 




54 


) 


489 


43 


RCL 


540 


01 


1 


54 


) 


490 


12 


12 


541 


54 


) 


54 


) 


491 


54 


) 


542 


33 




54 


) 


492 


33 


x'^ 


543 


53 


) 


42 


STO 


493 


54 


) 


544 


85 


+ 


08 


08 


494 


75 


- 


545 


53 


( 


99 


PRT 


495 


01 


1 


546 


43 


RCL 


91 


R/S 


496 


54 


) 


547 


08 


08 


76 


LBL 


497 


33 


x^ 


548 


33 


x2 


14 


D 


498 


54 


) 


549 


54 


) 


53 


( 


499 


85 


+ 


550 


54 


) 


53 


( 


500 


53 


( 


551 


54 


) 


04 


4 


501 


43 


RCL 


552 


42 


STO 


65 


X 


502 


08 


08 


553 


10 


10 


89 


IT 


503 


33 


x^ 


554 


99 


PRT 


65 


X 


504 


54 


) 


555 


66 


PAU 


43 


RCL 


505 


54 


) 


556 


53 


( 


11 


11 


506 


54 


) 


557 


43 


RCL 


54 


) 


507 


42 


STO 


558 


09 


09 


65 


X 


508 


09 


09 


559 


75 


- 


53 


( 


509 


99 


PRT 


560 


43 


RCL 


53 


( 


510 


91 


R/S 


561 


10 


10 


43 


RCL 


511 


76 


LBL 


562 


54 


) 


08 


08 


512 


15 


E 


563 


42 


STO 


65 


X 


513 


53 


( 


564 


11 


11 


01 


1 


514 


53 


( 


565 


99 


PRT 


05 


5 


515 


04 


4 


566 


91 


R/S 


00 


0 


516 


65 


X 








00 


0 


517 


89 


IT 








00 


0 


518 


65 


X 








00 


0 


519 


53 


( 








55 


T 


520 


43 


RCL 








53 


( 


521 


11 


11 








89 


IT 


522 


33 


x" 








65 


X 


523 


54 


) 








02 


2 


524 


54 


) 








65 


X 


525 


55 


T 








43 


RCL 


526 


53 


( 








12 


12 


527 


53 


( 








54 


) 


528 


53 


( 








54 


) 


529 


53 


( 








54 


) 


530 


53 


( 








55 




531 


43 


RCL 








53 


( 


532 


07 


07 








53 


( 


533 


55 


T 








53 


( 


534 


43 


RCL 








53 


( 


535 


12 


12 








53 


( 


536 


54 


) 
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APPENDIX D 



NUMERICAL INTEGRATION TI 59 PROGRAM 



A. INTRODUCTION 

In order to perforin the numerical integration of an (a) da, 
a standard Texas Instrument's program was used [Ref. 25]. 

This program performs the integration by using Simpson's 
discrete approximation based on the following expression 

X , 

I = / ^ f(x)dx ,(f +4f i+2f2+4f 3+2f4+ +2f _^+4f +f ) 

^ .j o n ii"“ 1 



where 



f(x) must be known at n+1 equally spaced points * 



h = 



X -X 

n o 
n 



n o 



n+7 = number of data registers available 

n = number of subintervals = 2, 4, 6, .... 



B. PROGRAM INSTRUCTIONS 

The program is taken from the master library program 
package by using the code 2nd Pgm 10 on the calculator. 
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The user's instructions are as follows 



Procedure 


Enter 


Press 


Display 


Enter data 


Subintervals 


A 


nl 


Enter data 


h 


B 


h^ 


Enter data 


Function values 








0 


C 


0 




fo 


R/S 


^o^ 




fl 

• 


R/S 


^1^ 




• 

• 

f 

n 


♦ 

R/S 


f ^ 
n 


Calculate 


— 


D 





^These values are printed automatically if the calculator 
is connected to PC-IOOA Print Cradle. 



The program steps are listed below giving location (LOG) , 



key symbol (KEY) , and comments. 



000 


76 


LBL 


019 


99 


PRT 


038 


06 


6 


001 


11 


A 


020 


92 


RTN 


039 


54 


) 


002 


53 


( 


021 


76 


LBL 


040 


42 


STO 


003 


50 


1x1 


022 


52 


EE 


041 


01 


01 


004 


42 


STO 


023 


00 


0 


042 


32 


x^t 


005 


05 


05 


024 


35 


1/x 


043 


98 


ADV 


006 


55 


T 


025 


92 


RTN 


044 


92 


RTN 


007 


02 


2 


026 


76 


LBL 


045 


76 


LBL 


008 


54 


) 


027 


12 


B 


046 


50 


1x1 


009 


42 


STO 


028 


42 


STO 


047 


76 


ST* 


010 


02 


02 


029 


03 


03 


048 


01 


01 


Oil 


22 


INV 


030 


99 


PRT 


049 


32 


x^t 


012 


59 


INT 


031 


92 


RTN 


050 


01 


1 


013 


29 


CP 


032 


76 


LBL 


051 


44 


SUM 


014 


22 


INV 


033 


13 


C 


052 


01 


01 


015 


67 


EQ 


034 


53 


( 


053 


32 


x>t 


016 


52 


EE 


035 


24 


CE 


054 


99 


PRT 


017 


43 


RCL 


036 


85 


+ 


055 


92 


RTN 


018 


05 


05 


037 


32 


x< t 


056 


61 


GTO 



122 



057 


50 


1x1 


108 


44 


SUM 


058 


76 


LBL 


109 


04 


04 


059 


14 


D 


110 


53 


( 


060 


53 


( 


111 


43 


RCL 


061 


43 


RCL 


112 


03 


03 


062 


05 


05 


113 


55 


T 


063 


85 


+ 


114 


03 


3 


064 


06 


6 


115 


54 


) 


065 


54 


) 


116 


49 


PRD 


066 


42 


STO 


117 


04 


04 


067 


01 


01 


118 


43 


RCL 


068 


73 


RC* 


119 


04 


04 


069 


01 


01 


120 


98 


ADV 


070 


42 


STO 


121 


99 


PRT 


071 


04 


04 


122 


92 


RTN 



072 76 LBL 

073 45 

074 01 1 

075 22 INV 

076 44 SUM 

077 01 01 

078 53 ( 

079 73 RC* 

080 01 01 

081 65 X 

082 04 4 

083 54 ) 

084 44 SUM 

085 04 04 

086 01 1 

087 22 INV 

088 44 SUM 

089 01 01 

090 22 INV 

091 97 DSZ 

092 02 02 

093 33 x^ 

094 53 ( 

095 73 RC* 

096 01 01 

097 65 X 

098 02 2 

099 54 ) 

100 44 SUM 

101 04 04 

102 61 GTO 

103 45 

104 76 LBL 

105 33 x^ 

106 73 RC* 

107 01 01 
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APPENDIX E 



TURN RATE LIMITATION TI 59 PROGRAMS 



A. INTRODUCTION 

The turn rate limitation calculations are divided into 
two parts with separate programs. 



B. PART I PROGRAM 

For a given initial detection range and a given ship 

speed. Part I performs a sequential calculation of 

<p„,K„,(p* and K* for each of the maximim turn rates investi- 
B B 

gated, together with the "no limit" conditions and 

based on 



^ . -1 ,3pld)imax 

4)^ = sin ' -90 



= 180 - - 4 )^ 



K 



A,B 



Ygd+cos 4 )^^b^ 
l^'max^^^^ ^A,B^ 



P-2 



= cos ^(P-2) 



V (l+P/2)^[l-(P/2)^]^"^/^ 



K* = 



4 > 



max 
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For p^2, the solution for * and K* are not valid. 
Furthermore, for solutions are valid for any 

of the quantities. 

The program is based on the fixed torpedo speed of 
= 35 kts (18 m/s) and an initial detection range of 

’'o = ' ' 250 m. 

The user instructions for the program are as follows: 



Prodcedure 


Enter 


Press 


Display 


Enter data 


Maximum Turn Rate 


A 


1 1 ^ 
' ^ ' max 


Enter data 


Ship speed 


B 


V 1 
T 


Calculate 
and K* 


— 


C 




Calculate cf). 
and 


— 


D 


A' ® 


Calculate K 
and Kg 


— 


E 





^These values are printed automatically if the calculator 
is connected to PC-IOOA Print Cradle. 



A listing of the program steps follows giving location 



(LOC) , code (COD) , key symbol (KEY) , and comments. 



000 


76 


LBL 


012 


69 


OP 


024 


54 


) 


001 


11 


A 


013 


00 


00 


025 


42 


STO 


002 


42 


STO 


014 


03 


3 


026 


03 


03 


003 


01 


01 


015 


03 


3 


027 


69 


OP 


004 


99 


PRT 


016 


69 


OP 


028 


06 


06 


005 


91 


R/S 


017 


04 


04 


029 


25 


CLR 


006 


76 


LBL 


018 


53 


( 


030 


69 


OP 


007 


12 


B 


019 


01 


1 


031 


00 


00 


008 


42 


STO 


020 


08 


8 


032 


02 


2 


009 


02 


02 


021 


55 


T 


033 


06 


6 


010 


99 


PRT 


022 


43 


RCL 


034 


03 


3 


Oil 


25 


CLR 


023 


02 


02 


035 


01 


1 
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036 


02 


2 


037 


07 


7 


038 


69 


OP 


039 


04 


04 


040 


53 


( 


041 


43 


RCL 


042 


03 


03 


043 


65 


X 


044 


03 


3 


045 


54 


) 


046 


42 


STO 


047 


10 


10 


048 


69 


OP 


049 


06 


06 


050 


25 


CLR 


051 


69 


OP 


052 


00 


00 


053 


04 


4 


054 


02 


2 


055 


03 


3 


056 


06 


6 


057 


03 


3 


058 


01 


1 


059 


02 


2 


060 


07 


7 


061 


69 


OP 


062 


04 


04 


063 


53 


( 


064 


53 


( 


065 


53 


( 


066 


03 


3 


067 


65 


X 


068 


01 


1 


069 


08 


8 


070 


54 


) 


071 


65 


X 


072 


43 


RCL 


073 


01 


01 


074 


54 


) 


075 


34 


/x 


076 


54 


) 


077 


42 


STO 


078 


11 


11 


079 


69 


OP 


080 


06 


06 


081 


91 


R/S 


082 


76 


L3L 


083 


13 


C 


084 


25 


CLR 


085 


69 


OP 


086 


00 


00 



03 


3 


138 


45 




03 


3 


139 


43 


RCL 


02 


2 


140 


03 


03 


03 


3 


141 


54 


) 


05 


5 


142 


54 


) 


01 


1 


143 


65 


X 


69 


OP 


144 


53 


( 


04 


04 


145 


53 


( 


53 


( 


146 


01 


1 


53 


( 


147 


75 


- 


53 


( 


148 


53 


( 


43 


RCL 


149 


53 


( 


03 


03 


150 


43 


RCL 


55 


T 


151 


03 


03 


02 


2 


152 


55 


T 


54 


) 


153 


02 


2 


22 


INV 


154 


54 


) 


39 


COS 


155 


33 


x" 


54 


) 


156 


54 


) 


54 


) 


157 


54 


) 


42 


STO 


158 


45 


yX 


04 


04 


159 


53 


( 


69 


OP 


160 


01 


1 


06 


06 


161 


75 


- 


25 


CLR 


162 


53 


( 


69 


OP 


163 


43 


RCL 


00 


00 


164 


03 


03 


02 


2 


165 


55 


f 


06 


6 


166 


02 


2 


05 


5 


167 


54 


) 


01 


1 


168 


54 


) 


69 


OP 


169 


54 


) 


04 


04 


170 


54 


) 


53 


( 


171 


55 


T 


53 


( 


172 


43 


RCL 


53 


( 


173 


01 


01 


53 


( 


174 


54 


) 


43 


RCL 


175 


42 


STO 


02 


02 


176 


05 


05 


65 


X 


177 


69 


OP 


53 


( 


178 


06 


06 


53 


( 


179 


91 


R/S 


01 


1 


180 


76 


LBL 


85 


+ 


181 


14 


D 


53 


( 


182 


25 


CLR 


43 


RCL 


183 


69 


OP 


03 


03 


184 


00 


00 


55 


4 


185 


03 


3 


02 


2 


186 


03 


3 


54 


) 


187 


02 


2 


54 


) 


188 


03 


3 



087 

088 

089 

090 

091 

092 

093 

094 

095 

096 

097 

098 

099 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 



126 



189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 



01 


1 


240 


06 


06 


291 


54 


) 


03 


3 


241 


91 


R/S 


292 


54 


) 


69 


OP 


242 


76 


LBL 


293 


54 


) 


04 


04 


243 


15 


E 


294 


42 


STO 


53 


( 


244 


25 


CLR 


295 


08 


08 


53 


( 


245 


69 


OP 


296 


69 


OP 


53 


( 


246 
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C. PART II PROGRAM 

Given K from Part I, Part II program performs the probe 
calculations of the corresponding initial approach angles 
^*^o^ A' ^*^0^ *' relationship 



The program used a fixed torpedo speed of = 35 kts (18 m/s) 



of this general probe program is to locate roots of the given 
function y = £(<P^) to evaluate the slope of the tangent line, 
and to find the maximum and minimum points on a graph. We 
will only use the first feature. 




P-1 



-r [ (1+ (cosd) 
o 




and an initial detection range of r^ = R^^ = 250 m. 



The purpose 
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The program requires a subroutine for the function to 
be investigated. This subroutine starts at program location 
140 and is located at label 2nd C. The user's instructions 
for the program are as follows : 



Procedure 


Enter 


Press 


Display 


Enter data 


Starting Value 
of (j) 

o 


2nd D 


^o 


Enter data 


The increment Acb 

o 


2nd E 


A4.0 


Calculate 

f(?^+A4.o> 


— 


A 


f(<^^+Acj,o) 


Calculate 


— 


B 


f (<|,o-A^o^ 


Display current 
value of 

o 


— 


E 


‘i’o 



If the value of Ac})^ is chosen too large, A<})^ may be 
replaced by A(})^/10 by pressing label 2nd A. Similarly, if 
a larger value of Atj)^ is required, Acj)^ can be replaced by 
10A<ii^ by pressing label 2nd B. A listing of the program 
steps follows, giving location (LOC) , code (COD) , key symbol 
(KEY), and comments. 

Associated with the subroutine, it should be noted that 
the value of K,r^ and p are entered separately in the 
memory locations 10, 11, and 12, respectively. 
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APPENDIX F 



FIGURES 
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Fig. 2. Location of Weather Stations Along the Norwegian 
Coast . 
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Fig. 3. Average Monthly Wind Speed in Beaufort from the Weather Station "Andenes. 
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Fig. 4. Histogram of Significant Wave in Percent per Year 
from the Weather Station "Andenes." 
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Fig. 5. Resonant Bubble Density at 12 kHz as a Function of 
Depth. 
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lant Bubble Density at 38 kHz as a Fxanction of 
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Fig. 7. Resonant Bubble Density at 120 kHz as a Function 
of Depth. 
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Fig. 8. Ambient Noise Level Curves. 
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Typical Sound Speed Profiles in Norwegian Coastal 
Waters. 
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Range (m) 
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Fig. 12. 



Frequencies of Wind Directions in Percent for 
Stations Along the Coast from "Hillesoy" to "Ona. " 
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Fig. 13. Frequencies of Wind Directions in Percent for 

Stations Along the Coast from "Myken" to "Furuholmen. " 
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Fig. 14. 
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Absorption Coefficient in dB/m as a Function of 
Frequency. 

145 



Source Level (dB re 1 uPa) 



SLi 



120 - 



no - 



too - 



9o 



Fig. 15. 
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Average Radiated Spectriam Level for Surface Ship 
as a Function of Speed in kts. 
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Fig. 16. General Noise Spectrum for a Cavitating Propeller. 
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Fig. 17. Spectrxun Level as a Fionction of Total Gas Content 
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Fig. 18. Measured Model and Full Scale Noise Spectra. 
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Fig. 20. Required Input S/N Ratio and BT Product for Various 
Operating Probabilities . 
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Fig. 21. Probability of Detection Versus Input S/N Ratio for 
Various Thresholds. 
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Fig. 22. -201ogR -aR Versus R for Frequency of 60 kHz. 



153 



Range (m 



30 kH:t 




Fig. 23. -201ogR -aR Versus R for a Frequency of 30 kHz. 
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Fig. 24a. Geometry at the Sea Surface Scattering. 
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Fig. 24b. Specular Scattering Geometry. 
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Fig. 25. 



Rough Surface Shadowing Geometry. 
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Fig. 26. Calculation Scheme of the Surface Scattering Effect. 
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Fig. 27. AIL Versus 



Fig. 




28. Geometry of the Sxibsurface Propagation Model. 
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Fig. 29. a) 



Chains of Successive Scattering. 



b) Scattering Pattern Going Through the Same 
Scatterer More than Once. 
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Fig. 30. Scattering Pattern. 
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Fig. 32. a and for 30 kHz. 
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Fig. 33. Interpolated Bubble Data. 
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Fig. 34. a^(a)n(a)da for 60 kHz. 
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Fig. 40. Hit Criterion Geometry. 
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Fig. 41. AOB-limitation Versus Ship Speed with Parameter 
Maximum Turn Rate. 
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